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ABSTRACT 

Using light curves from the HATNet survey for transiting extrasolar planets we investigate the 
optical broad-band photometric variability of a sample of 27, 560 field K and M dwarfs selected by 
color and proper-motion (V — K > 3.0, /x > 30 mas/yr, plus additional cuts in J — H vs. H — Kg 
and on the reduced proper motion) . We apply a variety of variability selection algorithms on the light 
curves to search for periodic and quasi-periodic variations, and for large-amplitude, long-duration 
flare events. To set the selection thresholds we conduct Monte Carlo simulations of light curves with 
realistic noise properties. A total of 3496 stars exhibit potential variability, including 95 stars with 
eclipses and 60 stars with flares. Based on a visual inspection of these light curves and an automated 
blending classification, we select 1928 stars, including 79 eclipsing binaries, as secure variable star 
detections that are not high probability blends. We find that only 43 of these stars, including 7 
of the eclipsing binaries, have previously been identified as variables or are blended with previously 
identified variables. One of the newly identified eclipsing binaries is 1RXS J154727. 5+450803, a known 
P = 3.55 day, late M-dwarf SB2 system, for which we derive preliminary estimates for the component 
masses and radii of M x = M 2 = 0.258 ± 0.008 M and Ri = R 2 = 0.289 ± 0.007 R . The radii of 
the component stars are larger than theoretical expectations if the system is older than ~ 200 Myr. 
The majority of the variables are heavily spotted BY Dra-type stars for which we determine rotation 
periods. This is the largest sample of photometric rotation periods for field K and M dwarfs published 
to date. Using this sample, we investigate the relations between period, color, age, and activity 
measures, including optical flaring, for K and M dwarfs. We find that the fraction of stars that are 
variable with peak-to-peak amplitudes greater than 0.01 mag increases exponentially with the V — Kg 
color such that approximately half of field dwarfs in the solar neighborhood with M < 0.2 Mq are 
variable at this level. Our data hints at a change in the rotation-activity-age connection for stars with 
M<0.25M Q . 

Subject headings: stars: late-type, low-mass, rotation, flare, fundamental parameters (masses, radii) 
— binaries: eclipsing — surveys — catalogs — techniques: photometric 



1. INTRODUCTION 

With the advent of large-scale automated time-series 
photometric surveys to search for microlensing events, 
transiting planets, and supernovae, well-sampled light 
curves for millions of stars have been, and are con- 
tinuing to be, collected. While the primary scientific 
goals of these surveys are to search for rare phenom- 
ena, the enormous databases of light curves that are a 
by-product of these surveys presents a great opportunity 
to study many other top ics related to stellar variabil- 
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data produced by dedicated wide-field transit surveys has 
been relatively under-utilized for studying topics unre- 
lated to transiting planets. In this paper we use data 
from the Hung arian-made Autom ated Telescope Net- 
work (HATNet; iBakos et all 12004 ) project to study the 
variability of probable K and M dwarf field stars. 

Combining photometric observations with proper mo- 
tion measurements is an effective method for selecting 
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nearby dwarf stars. This technique is routinely used in 
the search for cool st ars in the solar neighborhood (e.g. 
iRevle fc Ro bin 2003), and has been suggested as an ef- 
fective method for scree ning giants from transit surveys 
(| Gould fe Morgan! [20031 By selecting red, high proper 
motion stars, it is possible to obtain a sample that con- 
sists predominately of nearby K and M-dwarfs, with very 
few luminous, distant giants. Whereas a general vari- 
ability survey of Galactic field stars will yield a melange 
of objects that are o ften difficult to classify without de- 
tailed follow-up fe.g. lHartman et aL| [2004: Sh porer et al.1 
2007), by focusing a survey on red, high proper motion 
stars, one can narrow in on a few specific topics related 
to low-mass stars. 

Main-sequence stars smaller than the Sun are not 
known to exhibit significant pulsational instabilities, but 
they may exhibit the following types of photometric vari- 
ability: 1. Variability due to binarity (either eclipses or 
proximity effects such as ellipsoidal variability). 2. Vari- 
ability due to the rotational modulation or temporal evo- 
lution of starspots. 3. Flares. We discuss each type of 
variability, and what might be learned from studying it, 
below. 

1.1. Low- Mass Eclipsing Binaries 

In recent years there has been an increasing number 
of EBs discovered with K and M dwa rf components (see 
the list of 13 such binaries compiled bv lDevor et al.ll2008l . 
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note that most of these were found in the last 5 years). 
From these discoveries it has become clear that the radii 
of early M dwarfs and late K dwarfs are somewhat larger 
than p redicted by theory (the number typically stated 
is 10%;[lbrpr fe Ribas 2002; Ribas 2003; Lopez-Moralel 
|2005j [Ribas 2006; Bcat ty et al1l2007t lFernandedl2009f ) 

Most of the binaries found to date have periods shorter 
than a few days and are expected to have rotation peri- 
ods that are tidally synchronized to the orbital period. 
The rapid rotation in turn yields enhanced magnetic ac- 
tivity on these stars compared to isolated, slowly rotat- 
ing stars. It has been suggested that the discrepancy 
between theory and observation for these binary star 
components may be due to th eir enhanced magnetic ac- 
tivity inhibiting convection. (|Ribasl [2006t iTorres et al.l 
1200a lLopez Morales! [2007t iChabrier et all 120071) . Sui> 
port for this hypothesis comes from interfcromctric mea- 
surements of the luminosity-radius relation for inactive 
single K and M dwarfs which ap pears to be in agree - 
ment with theoretical predictions (jDemorv et al.l 12009). 
There is also evidence th at the discrepancy m ay be cor- 
related with metallicity (jLopez-MoralesI [20071 ) . Testing 
these hypotheses will require finding additional binaries 
spanning a range of parameters (mass, rotation/orbital 
period, metallicity, etc.). 

As transiting planets are discovered around smaller 
stars, the need for models that provide accurate masses 
and radii for these stars has become acute. For example, 
the errors in the planetary parameters for the transiting 
Super-Neptune GJ 436b appear to be dominated by the 
uncertainties in the stellar pa rameters of the 0.45 M M- 
dwarf host star (|Torresll2007l note that the author gives 
errors in the mass and radius for the star of ~ 3% as- 
suming that the theoretical models for the luminosity of 
M-dwarfs are accurate while making no such assumption 
about the radius). Having a large sample of low mass 
stars with measured masses and radii would enable the 
determination of precise empirical relations between the 
parameters for these stars. 

1.2. Stellar Rotation 

The rotation period is a fundamental measurable prop- 
erty of a star. For F, G, K and early M main sequence 
stars there is a well-established relat ion between rota- 
tion, magnetic activity and age (e.g. | Skumanichl 119721 : 
iNoves et al.lll984T iPizzolato et al.ll2003f) . In addition to 
illuminating aspects of stellar physics, this relation in 
practice provides the best meth od for measuring the ages 
of field main sequ ence stars ((Mamaick & Hillenbrand 
l2rMlBarneg[2007l) . 

For late M-dwarfs the picture is less clear. From a theo- 
retical standpoint one might expect that fully convective 
stars should not exhibit a rotation-activity connection if 
the connection is due to the aO-dynamo process, which 
operates at the interface of the radiative and convective 
zones in a star, and is thought to generate the large scale 
ma gnetic fields in the Sun ( see for example the discussion 
bv lMohantv fc Basril 120031 ). Nonetheless, several stud- 
ies have found evidence that the rotatio n-activity con- 
nection continues to very late M-dwarfs (iDelfosse et al.l 
H998tlMohanty & Basrill2003t iReiners fc Basrill2007h . In 
these studies the rotation period is inferred from the pro- 
jected rotation velocity vsini, which is measured spec- 
troscopically, while the degree of magnetic activity is es- 



timated by measuring either the Ha emission or X-ray 
emission. Rotation studies of this sort suffer both from 
the inclination axis ambiguity, and from low sensitivity 
to slow rotation. In practice it is very difficult to mea- 
sure vsmi values less than ~ 1 km/s, which generally 
means that it is only possible to place lower limits on 
the period for late M-dwarf stars with periods longer 
than ~ 10 days. Moreover, because low-mass stars are 
intrinsically faint, these studies require large telescopes 
to obtain high-resolution, high S/N spectra, so that typ- 
ically only a few tens of stars are studied at a time. 

There are two techniques that have been used to di- 
rectly measure stell ar rotation periods. Th e first tech- 
nique, pioneered bv lWilson &; Bappul (|1957l ). is to mon- 
itor the emission from the cores of the Call H and K 
lines, searching for periodic variations. The venerable 
Mount-Wilson Observatory HK project has used this 
technique to measure the rotation periods of more than 
100 slowly rotating dwarfs and giant stars (lWilsonlll978l : 
iDuncan et all 119911 : iBaliunas et aH 1 19951 [1996J). Alter- 
natively, if a star has significant spot-coverage it may 
be possible to measure its rotation period by detecting 
quasi-periodic variations in its broad-band photometric 
brightness. Studies of this sort ha ve been carried out 
in abundance for open clusters (e.g. iRadick et al.lll987l) 
as w ell as for some field stars (e.g. IStrassmeier et all 
2000). While there are rich samples of rotation pe- 
riods for K and M dwarfs in open clusters with ages 
< 600 Mvr (llrwin et all 2006. 2007L l2009aHMeibom et all 
120091 : lHartman et al.ll2009f ). the data for older K and M 
dwarfs is quite sparse. As such, there are few obser- 
vational constraints on the rotational evolution of these 
stars after ~ 0.5 Gyr. 

Unlike spectroscopic studies, photometric surveys may 
yield rotation periods for hundreds of stars at a time. 
There are, however, some drawbacks to these surveys. 
The spot distribution on the surface of a star may in 
general be quite complex, so the resulting signal in the 
light curve will not always take a simple form. Since the 
number of brightness minima per cycle is not known a 
priori, there is a risk that the true rotation period may 
be a harmonic of the measured period. Spots on the 
Sun come and go on time-scales shorter than the Solar 
rotation period, and indeed stellar light curves also ex- 
hibit secular trends. For long period stars the measured 
variation time-scale may actually correspond to a spot 
evolution time-scale rather than the rotation period of 
the star. For short period stars there may be difficulties 
in distinguishing spot modulation from binarity effects 
(though for these stars the rotation period is expected to 
be tidally synchronized to the orbital period). 

Despite these caveats, given the existing uncertainties 
in the rotation-activity connection for low mass stars and 
the potential to use rotation as a proxy for age, a large, 
homogeneously collected sample of photometric rotation 
periods for field K and M dwarfs could potentially be of 
high value. 

1.3. Flares 

Flaring is known to be a common phenomenon among 
K and M dwarfs. Studies of open cluster and field 
flare stars have shown that the frequency of flaring 
increases with decreasing stellar mass, and decreases 
with increasing stellar age (e.g. lAmbartsumvan et al.l 
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Il970t iMirzovan et al.lll989T) . Significant flaring on these 
low-mass dwarfs is likely to impac t the habitabi l ity of 
any planets they may harbor (e.g. iKasting et alj [19931 ; 
iLammer et alj2007t iGuinan fc Engld2009[ ) sodetermin- 
ing the frequency of flares, and its connection with other 
stellar properties such as rotation, has important impli- 
cations for the study of exoplanets. 

1.4. The HATNet Survey 

To address these topics we use data from HATNet to 
conduct a variability survey of K and M dwarfs. The on- 
going HATNet project is a wide-field search for transit- 
ing extrasolar planets (TEPs) orbiting relatively bright 
stars. The project employs a network of 7 robotic tele- 
scopes (4 in Arizona at Fred Lawrence Whipple Obser- 
vatory, 2 in Hawaii on the roof of the Sub-Millimeter 
Array at Mauna Kea Observatory, and 1 in Israel at 
Wi se Observatory; the latter is referred to as WHAT, 
see lShporer et all 12009 ) which have been used to obtain 
some ~ 700,000 images covering approximately 10% of 
the sky. The survey has generated light curves for ap- 
proximately 2.5 million stars, from which ~ 900 candi- 
date TEPs have been identified. To date, the survey has 
announ ced the discovery of 12 TEPs, including HAT- 
P-llb (jBakos et al J 12009 ), a Super-Neptune (0.08 Mj) 
planet that is the smallest found so far by a ground-based 
transit survey. While the primary focus of the HATNet 
project has been the discovery of TEPs, some results not 
related to planets have also been presented. This includes 
the discovery and analysis of a low-mass M dwarf in a 
single -lined eclipsing binary (EB) system (|Beattv et al.l 
I2007D , searches for variable stars in two HATNet /WHAT 
fields (|Hartman et al.ll2004t IShporer et al.ll2007f ). 

1.5. Overview of the Paper 

The structure of the paper is as follows. In § [2] we 
describe both the HATNet photometric data, and select 
the sample of field K and M dwarfs. In § [3] we discuss 
our methods for selecting variable stars. We estimate 
the degree of blending for potential variables in § H) We 
match our catalog of variables to other catalogs in § [5j 
We discuss the properties of the variables in §|5] including 
an analysis of one of the EB systems found in the survey. 
We conclude in §[71 In appendix lAlwe describe the Monte 
Carlo simulations used in establishing our variability se- 
lection thresholds, while in appendix [B] we present the 
catalog of variable stars. 

2. OBSERVATIONAL DATA 

2.1. HATNet Data 

The HATNet project, which has been in operation 
since 2003, uses a network of 7 small (11cm aperture), 
autonomous telescopes to obtain time-series observations 
of star s. For details on the system design see lBakos et al.1 
( 2004) ; here we briefly review a few points that are rele- 
vant to the survey presented in this paper. Prior to 2007 
September each t elescope emplo yed a 2K x 2K CCD and a 
Cousins I c filter (|Cousinslll97fil ). The 2Kx2K CCDs cov- 
ered an 8.2° x 8.2° field of view (FOV) at a pixel scale of 
14". With these CCDs, stars with 7.5 < I c < 14.0 were 
observed with a typical per-image photometric precision 
of a few mmag at the bright end, 0.01 mag at Ic ~ 11, 
and 0.1 mag at Ic ~ 13.5. After this date the telescopes 



were refitted with 4Kx4K CCDs and Cousins Rc filters. 
The new CCDs cover a 10.6° x 10.6° FOV at a pixel scale 
of 9". With these CCDs, stars with 8.0 < R c < 15.0 are 
observed with a typical per-image photometric precision 
of a few mmag at the bright end, 0.01 mag at Rc ~ 12, 
and 0.1 mag at Rc ~ 15. The exact magnitude lim- 
its and precision as a function of magnitude vary within 
a field due to vignetting, and from field to field due to 
differences in the reduction procedure used and in the 
degree of stellar crowding. In 2008 September the filters 
were changed to Sloan r, though we do not include any 
observations taken through the new filters in the survey 
presented here. For both CCD formats, the typical full 
width at half maximum (FWHM) of the point spread 
function (PSF) is ~ 2 pixels (i.e. - 30" for the 2K fields 
and - 20" for the 4K fields). 

The data for this survey comes from 72 HATNet fields 
with declinations between +15° and +52°. These fields 
are defined by dividing the sky into 838 7.5° x 7.5° tiles. 
The survey covers approximately 4000 square degrees, or 
roughly 10% of the sky. 

The data reduction pipeline has evolved over time, as 
such the fields studied in this survey have not all been 
reduced in a uniform manner. For simplicity we choose 
to use the available light curves as is, rather than re- 
reducing the fields in a consistent manner optimized for 
finding variable stars rather than TEPs. Both aper- 
ture photometry (AP) and image subtraction photome- 
try (ISM) have been used for reductions. Both pi pelines 
were developed from scratch for HATNet. SeelPall (|2009 ) 
for detailed descriptions of both methods. 

For bo th pipelines the Tw o Micron All-Sky Survey 
(2MASS: ISkrutskie et al.ll2006l ) is used as the astrometric 
reference. The astrometric solutions for t he images are 
determ ined using the methods described bv lPal fc Bakos 
(|2006t ) and EH (|2009 ). Photometry is performed at the 
positions of 2MASS sources transformed to the image 
coordinate system. For each resulting light curve the 
median magnitude is fixed to the Ic or Rc magnitude of 
the source based on a transformation from the 2MASS 
J,H and K$ magnitudes. For the ISM reduction this 
magnitude is also used as the reference magnitude for 
each source in converting differential flux measurements 
into magnitudes. 

Both the AP and ISM pipelines produce light curves 
that are calibrated against ensemble variations in the flux 
scale (for AP this is done as a step in the pipeline, for 
ISM this is an automatic result of the method) . For each 
source, light curves are obtained using three separate 
apertures. The set of apertures used has changed over 
time; the most recent reductions use aperture radii of 
1.45, 1.95 and 2.35 pixels. Following the post-processing 
routines discussed below, we adopt a single "best" aper- 
ture for each light curve. 

The calibrated light curves for each aperture are passed 
through two routines that remove systematic variations 
from the light curves that are not corrected in cali- 
brating the ensemble. The first routine (EPD) decor- 
relates each light curve against a set of external pa- 
rameters including parameters describing the shape of 
the PSF, the sub-pixel position of the star on the im- 
age, the zenith angle, the hour angle, the local sky 
background, and the variance of the background (see 
iBakos et ai1l2009) . After applying EPD, the light curves 
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are th en processed with the Trend- Filtering Algorithm 
(TFA; iKovacs. Bakos &: Novesl 120051 ) which decorrelates 
each light curve against a representative sample of other 
light curves from the field. The number of template light 
curves used differs between the fields, typically the num- 
ber is ~ 8% of the total number of images for that field. 
In applying the TFA routine we also perform er-clipping 
on the light curves since this generally reduces the num- 
ber of false alarms when searching for transits. For the 
remainder of the paper we will refer to light curves that 
have been processed through EPD only, without applica- 
tion of TFA, as EPD light curves, and will refer to light 
curves that have been processed through both EPD and 
TFA as TFA light curves. We note that for some fields 
the EPD light curves were not stored and only TFA light 
curves are available. 

Both of these algorithms tend to improve the signal 
to noise ratio of transit signals in the light curves, but 
they may distort the light curves of stars that show 
large-amplitude, long-period, continuous variability. Ad- 
ditionally the decorrelation against the zenith and hour 
angles in the EPD routine will tend to filter out real 
variable star signals with periods very close to a side- 
real day or an integer multiple of a sidereal day. The 
TFA routine in particular may distort long-period sig- 
nals while increasing the signal to noise ratio of short- 
period signals. For this reason we analyze both the EPD 
and TFA light curves, when available, to select variable 
stars (§ [3]). We note that for the analysis in this pa- 
per we do n ot use the signal-reconstructio n mode TFA 
presented bv lKovacs, Bakos fc Novesl (|2005|) . Once a sig- 
nal is detected, TFA can be run in this mode to obtain a 
trend-filtered light curve that is free of signal distortions, 
however for signal detection one must use general TFA 
since the signal is not known a priori. 

Finally, an optimal aperture is chosen for each star. 
For stars fainter than a fixed limit the smallest aperture 
is used (to minimize the sky noise), for brighter stars 
the aperture with the smallest root-mean-square (RMS) 
light curve is used. 

2.2. Composite Light Curves 

Because the separation between the HATNet field cen- 
ters is smaller than the FOV of the HATNet telescopes 
for both the 2K and 4K CCDs, some stars are observed in 
multiple fields. These stars may have more than one light 
curve, which we combine into composite light curves. In 
making a composite light curve we subtract the median 
magnitude from each component light curve. For fields 
reduced with both ISM and AP we use the light curve 
with the lowest RMS, and in the case of equal RMS we 
use the ISM light curve. Note that the composite light 
curve for a star may include a mix of 1q and Rc photom- 
etry. While the amplitude of variability may be different 
from filter to filter, the period and phasing for variations 
due to eclipses or the rotational modulation of starspots 
will be independent of bandpass. For simplicity we do 
not allow for independent amplitudes of different filters 
in searching for variability using the methods described 
in § [3] we do not expect this to make a significant dif- 
ference to period detections. However, we note that a 
simple merging of photometry from different filters may 
result in spurious side lobes in the power spectrum; for 
a more detailed analysis of individual objects this effect 



should be considered. 

2.3. Selection of the K and M Dwarf Sample 

To select the sample of stars that are probable K and 
M dwarfs we apply cuts on the proper motion and on 
the color. Proper motion measureme nts are taken from 
the PPM-Extended catalog fPPMX; Hoser et all [2001) 
which provides proper motions with precisions ranging 
from 2 mas/yr to 10 mas/yr for 18 million stars over the 
full sky down to a limiting magnitude of V ~ 15.2. The 
PPMX catalog provides complete coverage of the HAT- 
Net survey for stars with V — Ic ^$ 1-2 for the 2K fields 
and for stars with V - R < 0.7 for the 4K fields. For 
stars redder than these limits, the faint limit of HATNet 
is deeper than the faint limit of PPMX. We select all stars 
from this catalog with a proper motion [i > 30 mas/yr. 
For the color selection we use the 2MASS JHK$ pho- 
tometry and, where available, y-band photometry from 
the P PMX catalog, w hich is taken from the Tycho-2 cat- 
alog (iHpg et al.ll2000D a nd tr ansformed to the Johnson 
system by iKharchenkol (|2001l ). Only ~ 4% of the stars 
have V photometry given in the PPMX catalog, for the 
majority of stars that do not have V photometry we cal- 
culate an approximate V magnitude using 

V = -0.0053 + 3.5326 J + 1.314m - 3.8331X S . (1) 

which is deter mined from 590 Landolt Standard stars 
(jLandoldl 19921 ) with 2MASS photometry, and is used in- 
ternally by the HATNet project to estimate the V mag- 
nitudes of transit candidates for follow-up observations. 
We then select stars with V — K$ > 3.0 which corre- 
sponds roughly to star s with spectral types later than K6 
(jBessell fc Brett! 119881) . This selects a total of 471,970 
stars from the PPMX catalog, of which 33, 177 fall in a 
reduced HATNet field; of these 32,831 have a HATNet 
light curve containing more than 1000 points. 

Note that extrapolating a V magnitude from near- 
infrared photometry is not generally reliable to more 
than a few tenths of a magnitude. We therefore also 
obtained V magnitudes for star s in our sample by m atch- 
ing to the USNO-B1.0 catalog (jMonet et al.ll2003h . This 
matching was done after the variability search described 
in § [3J we choose not to redo the sample selection and 
the subsequent variability selection. Note that low-mass 
sub-dwarfs, which have anomalously blue V — K$ values, 
will pass a selection on V — K$ > 3.0 computed using 
cq. [T] to extrapolate the V magnitude, while they would 
not necessarily pass a selection using the measured value 
of V — Ks- To transform from the photographic Bjj, Rjj 
magnitudes in the USNO-B1.0 catalog to the 1^-band we 
use a relation of the form: 

V = aB v + bBf, + cR v + dR 2 v + eB v R v + f (2) 

with coefficients given separately in table [T] for the 

(Bu,i,Ru,i), (Bu,i,Ru,2), (Bu.2, Ru.i) and (Bjj,2,Bu,2) 
combinations. These transformations were determined 
using ~ 1100 stars with both V photometry in the 
PPMX catalog and USNO-B1.0 {B v> Ro) photometry. 
Based on the root-mean-square (RMS) scatter of the 
post-transformation residuals, we used the (-Bj/,2) Ru,2), 
(Bu,i,Ru,2), (Bu,2,Ru,i) and (Bu,i,Ru,i) transforma- 
tions in order of preference. For stars with neither 
PPMX V photometry nor USNO-B1.0 photometry, we 
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used eq. [TJ For the remainder of the analysis in this pa- 
per, the V magnitude is taken from PPMX (Tycho-2) 
for 3.1% of the stars in our sample, from USNO-B1.0 for 
93.6% of the stars, and is transformed from the 2MASS 
magnitudes for 3.3% of the stars. 

Figure [T] shows the J — H vs. H — K$ color-color dia- 
gram for the selected sample. We also show the expected 
relations for dwarf stars and for giants. T he relation for 
dwarf s is taken from a combination of the iBaraffe et al.l 
(1998) 1.0 Gyr isochrone for solar metallicity sta r s with 
0.15 M Q < M < 0.7 Mq and the IChabrier et all (|2000D 
models for objects with M < 0.075 M Q . The JHK mag- 
nitudes for t hese isochrones were tr ansformed from the 
CIT system (jElias et al.l 1 1981 fl98l to the 2 MASS sys- 
tem u sing the transformations determined by ICarpenterl 
(2001). The relation for giant stars with \ogg < 2.0 
is taken f rom the 1.0 Gyr, solar metallicity Padova 
isoch rone (|Marigo et al.l 120081 : iBonatto. Bica fe Girardil 
2004) 5 . While the majority of stars lie in the expected 
dwarf range, a significant number of stars fall along the 
giant branch. Some of these stars may be rare carbon 
dwarfs, but the majority are most likely giants with in- 
accurate proper motion measurements in the PPMX cat- 
alog. Of the 2445 selected stars with J - H > 0.8 that 
have HATNet light curves, 87% have undetected proper 
motions or proper motions less than 10 mas/yr in the 
USNO-B1.0 catalog, this is compared to 28% of the sam- 
ple with J — H < 0.8. A visual inspection of the POSS-I 
and POSS-II Digitized Sky Survey images for a number 
of the sources with J — H > 0.8 and /i > 100 mas/yr re- 
vealed none with visually detectable proper motion, and 
in many cases the object consists of two close, compa- 
rably bright stars, for which misidentification of sources 
may be to blame for the spurious proper motion detec- 
tion. This includes several stars where the PPMX and 
USNO-B1.0 proper motion values are comparable. We 
therefore apply an additional cut in the J — H vs H — K$ 
color-color diagram as shown in figure [1] to reduce the 
sample to 28, 785 stars. 

Figure [2] shows a V — J vs. Hj reduced proper- motion 
(RPM: [Luytenl[l922h diagram for the sample. Here the 
RPM, Hj, is calculated as 
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Hj = J + 51og 10 (/i/1000) 



(3) 



and gives a rough measure of the absolute magnitude Mj 
of a star. We show roughly the lines separating main 
sequence dwarfs from sub-dwarfs and giants. In figure [3] 
we compare the RPM to Mj for 239 stars in the sampl e 
which have a Hipparcos parallax (jPerrvman et al.lll997f ). 
To remove additional giants from the sample we reject 
1225 stars with Hj < 3.0, leaving our final sample of 
27, 560 stars. 

3. SELECTION OF VARIABLE STARS 

We use a number of techniques to search for 
light curves that show significant variability. 
The techniques include the phase-binning and 
harmonic-fitting Analysis of Variance periodqgram s 
(AoV/AoVHarm; ISchwarzenberg-Czernvl Il989l . |1996[) . 
the Discrete Auto-C orrelation Function (DACF; 
lEdelson fc Krolikl 1 1988T ) . and the Box-Least-Squares 

5 The isochrone was obtained from the CMD 2.1 web interface 
http: //stev. oapd. inaf . it /cgi-bin / cmd_2. 1 
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Fig. 1. — J — H vs. H — Kg color-color diagram for 32,831 
stars that have V — Kg > 3.0 with V taken either from the PPMX 
catalog (Tycho-2) or extrapolated from the 2MASS J, H and Kg 
magnitudes using eq. [l] fi > 30 mas/yr from the PPMX catalog, 
and that have a HATNet light curve containing more than 1000 
points (gray-scale poin ts). The solid line s hows the expected re - 
lation for cool dwarfs (Baraffc et al. 1998: Chabricr et al. 2000), 
while the dot-dashed line shows the expected re lation for giants 
UMarigo et al . 2008; Bonatt o~Bica fc G irardi 2003). Stars outside 
the area enclosed by the dotted line are rejected. For display pur- 
poses we have added slight Gaussian noise to the observed colors 
in the plot. 



E 




1 2 3 4 5 6 

V-J [mag] 

Fig. 2. — V — J vs. Hj RPM diagram for stars with > 
30 mas/yr passing V — K$ and JHK color cuts. Note that in this 
plot we use V magnitudes that are transformed from the USNO- 
B1.0 photographic magnitudes for the majority stars, and not the 
V magnitudes transformed from JHK that were used in applying 
the initial V — Kg cut. The lines separate main sequence dwarfs 
from sub-dwarfs and giants. We reject the 1225 stars with Hj < 
3.0. The lines separating main s equen ce dwarfs from sub-dwarfs 
are taken fromfYong & Lambert 1120031 ). 



(BLS; iKovacs. Zucker fc Mazehl l2002) alg orithms as im- 
plem ented in the VARTOOLS program 6 (Har tman et al.l 
2008). We also conduct a search for flare-like events in 
the light curves. We apply the algorithms to both the 
EPD and TFA light curves of sources, when available. 
For the flare-search we use only the EPD light curves 
because the cr-clipping applied to the TFA light curves 
may remove real flares. 

Because the light curves contain non-Gaussian, 
temporally-correlated noise, formal estimates for the 

6 The VARTOOLS program is available at 
http://www.cfa.harvard.edu/~jhartman/vartools/ 
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TABLE 1 

Coefficients for transformations from USNO-B1.0 Bu and Rjj magnitudes to V (eq.EJ. 



Bu 


Ru 


« 


6 


c 


d 


e 


/ 


RMS [mag] 


1 


1 


-0.76 ± 0.19 


-0.02 ± 0.01 


1.71 ± 0.15 


-0.14 ±0.01 


0.15 ±0.02 


1.63 ±0.52 


0.27 


1 


2 


0.72 ± 0.08 


-0.061 ±0.002 


0.38 ±0.06 


-0.057 ± 0.006 


0.114 ±0.007 


-0.89 ±0.26 


0.19 


2 


1 


0.85 ± 0.09 


-0.045 ± 0.006 


0.44 ± 0.08 


-0.051 ± 0.003 


0.079 ± 0.008 


-1.52 ±0.24 


0.22 


2 


2 


0.73 ±0.14 


-0.164 ±0.008 


0.46 ± 0.13 


-0.20 ±0.01 


0.35 ±0.02 


-0.75 ±0.26 


0.19 



at 

CO 




-5 5 
Mj [mag] 



10 



Fig. 3.— The absolute magnitude Mj vs. the RPM for 239 stars 
in the sample which have a Hipparcos parallax. This confirms that 
the RPM provides a rough estimate of the absolute magnitude for 
this sample of stars. We reject stars with Hj < 3.0, as these appear 
to be predominately giants. 

variability detection significance are unreliable. We 
therefore have conducted Monte Carlo simulations of 
light curves with realistic noise properties to inform our 
choice of selection thresholds for several of the algorithms 
mentioned above, these are described in appendix [Al In 
the following subsections we discuss the use of each of the 
variability selection algorithms in turn. We finish with 
a comparison of the different methods. The resulting 
catalog of variable stars is presented in appendix [B] 

3.1. AoV 

The AoV periodogram is a method for de- 
tecting continuous pe r iodic variations suggested by 
ISchwarzenberg-Czernvl (|1989h which typically yields 
higher S/N detections than oth er periodograms. The 
origin al method suggested by ISchwarzenberg-Czernvl 
(1989) uses phase-binning for the model signal, so it is 
most comparable to the p opular Phase Dispe rsion Mini- 
miz ation technique (PPM; Stcllingwerj ll978l ). Following 
this lSchwarzenberg-Czeriiv (|1996h introduced an efficient 
method for fitting a Fourier series to a non-uniformly 
sampled light curve using a basis of orthogonal poly- 
nomials. When combined with an ANalysis Of VAri- 
ance (ANOVA) statistic, the result is the AoVHarm pe- 
riodogram. We applied both methods to our light curves, 
and discuss each in turn. 

3.1.1. Phase Binning AoV - Search for General Periodic 

Variability 

We run the AoV algorithm with phase-binning on the 
full sample of stars. We first apply a 5er iterative clipping 
to the light curve before searching for periods between 
0.1 and 100 days. We use 8 phase bins, and generate the 
periodogram at a frequency resolution of 0.1/T where 



T is the total time-span covered by a given light curve. 
We then determine the peak at 10 times higher resolu- 
tion. As our figure of merit we use the signal-to-noise 
ratio (S/N), with a 5cr iterative clipping applied to the 
periodogram in calculating the noise (the RMS of the pe- 
riodogram). Figure 2] shows the AoV S/N as a function 
of the peak period for various light curve subsamplcs. 
For comparison we also show the Monte Carlo simulation 
results for each subset. We adopt a separate selection 
threshold on S/N for each subsample. The thresholds 
have the form 



S/N n 



f S/N if P < P 

{ S/N Q (P/P ) a if P>P 



(4) 



The adopted values of S/No, Pq and a are listed for each 
subsample in Tabled In addition to this selection we also 
reject detections with periods near 1 sidereal day, one of 
its harmonics, or other periods which appear as spikes 
in the histogram of detected periods (the latter includes 
periods between 5.71 and 5.80 days). For composite light 
curves that contain both 2K and 4K observations we take 

S/iVmin = f2I<S/N 1Jlin ^K + faC S / N m in,4K ( 5 ) 

where $ik and fiK are the fraction of points in the 
light curve that come from 2K and 4K observations, and 
S/N m i U; 2K and SyiVmin,4K are the 2K and 4K thresholds 
at the period of the composite light curve. 

As summarized in table [3[ our selection threshold 
passes a total of 1320 EPD light curves, and 1729 TFA 
light curves. These are inspected by eye to reject ob- 
vious false alarms and to identify EBs. There are 753 
EPD light curves that we judge to show clear, continu- 
ous, periodic variability, 47 that show eclipses, 419 that 
we consider to be questionable (these are included in the 
catalog, but flagged as questionable), and 101 that we 
reject. For the TFA light curves the respective num- 
bers are 1210, 64, 400, and 55. Note that the distinction 
between "clear" variability and "questionable" cases is 
fairly subjective. Generally we require the variations to 
be obvious to the eye for periods > 30 days, for shorter 
periods we consider the selection to be questionable if 
it appears that the variability selection may be due to 
enhanced scatter on a few nights. 

3.1.2. Harmonic AoV - Search for Sinusoidal Periodic 
Variability 

The AoVHarm periodogram is generated for each star 
in a similar manner to the AoV periodogram. We run the 
algorithm using a sinusoid model with no higher harmon- 
ics (it is thus comparable to D FT methods, or to the pop - 
ular Lomb-Scargle technique; Lomb 1976; Scargle 1982h . 
As for the phase-binning AoV search we use eq. [5] for the 
selection threshold, with parameters given for each sub- 
sample in table [2j Figure [5] shows the AoVHarm S/N as 
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TABLE 2 

Selection Threshold Parameters for AoV and AoVHarm 



AoV AoVHarm 



Chip 


LC Type 


Mag. 


S/N 


Po [days] 


Q 


S/N 


Po [days] 


a 


EPD 


2K 


6.5 


30.0 


5 


0.40 


40.0 


4 


0.41 


EPD 


2K 


7.5 


30.0 


5 


0.40 


40.0 


4 


0.41 


EPD 


2K 


8.5 


30.0 


5 


0.40 


40.0 


4 


0.41 


EPD 


2K 


9.5 


30.0 


5 


0.33 


30.0 


4 


0.37 


EPD 


2K 


10.5 


25.0 


5 


0.23 


30.0 


4 


0.37 


EPD 


2K 


11.5 


20.0 


5 


0.23 


25.0 


4 


0.43 


EPD 


2K 


12.5 


20.0 


5 


0.23 


25.0 


4 


0.36 


EPD 


2K 


13.5 


20.0 


5 


0.14 


25.0 


4 


0.22 


TFA 


2K 


6.5 


20.0 


5 


0.14 


30.0 


4 


0.16 


TFA 


2K 


7.5 


20.0 


5 


0.14 


30.0 


4 


0.16 


TEA 


2K 


8.5 


20.0 


5 


0.14 


30.0 


4 


0.16 


TFA 


2K 


9.5 


20.0 


5 


0.00 


20.0 


4 


0.22 


TFA 


2K 


10.5 


15.0 


5 


0.00 


20.0 


4 


0.13 


TFA 


2K 


11.5 


15.0 


5 


0.00 


20.0 


4 


0.00 


TFA 


2K 


12.5 


15.0 


5 


0.00 


20.0 


4 


0.00 


TFA 


2K 


13.5 


15.0 


5 


0.00 


20.0 


4 


0.00 


EPD 


4K 


7.5 


30.0 


5 


0.33 


30.0 


4 


0.50 


EPD 


4K 


8.5 


30.0 


5 


0.33 


30.0 


4 


0.50 


EPD 


IK 


9.5 


30.0 


5 


0.33 


30.0 


4 


0.50 


EPD 


4K 


10.5 


30.0 


5 


0.33 


30.0 


4 


0.50 


EPD 


4K 


11.5 


25.0 


5 


0.23 


25.0 


4 


0.43 


EPD 


4K 


12.5 


20.0 


5 


0.23 


25.0 


4 


0.36 


EPD 


4K 


12.5 


20.0 


5 


0.23 


20.0 


4 


0.43 


EPD 


4K 


13.5 


20.0 


5 


0.14 


20.0 


4 


0.28 


EPD 


4K 


14.5 


20.0 


5 


0.00 


20.0 


4 


0.22 


TFA 


4K 


7.5 


20.0 


5 


0.00 


25.0 


4 


0.15 


TFA 


4K 


8.5 


20.0 


5 


0.00 


25.0 


4 


0.15 


TFA 


4K 


9.5 


20.0 


5 


0.00 


25.0 


4 


0.15 


TFA 


4K 


10.5 


20.0 


5 


0.00 


25.0 


4 


0.06 


TFA 


4K 


11.5 


15.0 


5 


0.00 


20.0 


4 


0.00 


TFA 


4K 


12.5 


15.0 


5 


0.00 


20.0 


4 


0.00 


TFA 


4K 


12.5 


15.0 


5 


0.00 


20.0 


4 


0.00 


TFA 


4K 


13.5 


15.0 


5 


0.00 


20.0 


4 


0.00 


TFA 


IK 


14.5 


15.0 


5 


0.00 


20.0 


4 


0.00 



a function of peak period for several light curve subsam- 
ples. 

Our selection threshold passes a total of 1337 EPD 
light curves and 1717 TFA light curves. Inspecting these 
by eye we find 1082 EPD light curves that show clear, 
continuous, periodic variability, 34 that show eclipses, 
185 that we consider to be questionable, and 36 that we 
reject. For the TFA light curves the numbers are 1443, 
43, 217, and 14 respectively (see also tabled]). 

3.2. DACF - Search for Quasi- Periodic Variability 

The DACF is a technique that has frequently been used 
to determine the variation time-scale of quasi-periodic 
signals. In particular this technique is commonly ap- 
plied in measuring the rotation period of a star from 
a light curve that exhibits variations due to the rota- 
tional modulation of starspots which may be varying 
in size, intensity, p osition, or number over time (e.g. 
lAigrain et al.1 [2008). It is not obvious, however, that 
the DACF method provides better period determinations 
than Fourier methods such as AoVHarm even in quasi- 
periodic cases. Since even in cases where the coherence 
time-scale is short compared to the period, the Fourier 
power spectrum should still have a predominant peak 
near the period of the star that can be determined in a 
straightforward fashion, whereas the automatic identifi- 
cation of the predominant variation time-scale from an 
autocorrelation function is non-trivial, as seen below. It 
is, therefore, with some skepticism that we attempt to 
use the DACF method to identify periods. 



For each light curve we calculate the DACF at time lags 
ranging from to 100 days with a step-size of 1 day. The 
binning time-scale of 1 day means that we will only be 
sensitive to variation time-scales > 2 days. In practice, 
our peak finding algorithm limits us to periods > 10 days. 
A light curve with a significant periodic or quasi-periodic 
signal with a timescale T will have a DACF that peaks at 
T, as well as at 2T, 3T, . . . depending on the coherence 
of the signal. 

To automate the selection of peaks in the DACF we 
use the following routine (below, t/j is the DACF value 
for time-lag ti): 

1. Identify connected sets of points (ti, y{) with y$ > 0. 

2. Extend the left and right boundaries of each set 
until a local minimum is found in both directions. 

3. Reject sets with 3 or fewer points, with for the 
left boundary, or with the maximum time-lag com- 
puted for the right boundary. This leaves N pca k 
sets of points (peaks) to consider. 

4. Fit a quadratic function y(t) to each of the iV pea k 
sets. 

5. Letting Xn-3 an( i Xjv-i be the x 2 values from fit- 
ting a quadratic function and a constant function 
respectively to the set, we perform an F-test on the 
statistic 

r _ (Xn-i ~ Ajv-3)/^ / fi \ 
/_ X 2 N s/(N-3) W 
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TABLE 3 

Summary of variable star selection 



Variability 


LC Reduction 


Automatically 


Non-EB 








Selection Type 


Type 


Selected 


Var. 


EB 


Quest. 


Rej. 


AoV 


EPD 


1320 


753 


47 


419 


101 


AoV 


TFA 


1729 


1210 


64 


400 


55 


AoV 


Both a 


781 


576 


11 


53 


5 


AoV 


Either b 


2268 


1387 


70 


766 


151 


AoV only c 


Either b 


450 


120 


4 


317 


79 


AovHarm 


EPD 


1337 


1082 


34 


185 


36 


AovHarm 


TFA 


1717 


1443 


43 


217 


11 


AovHarm 


Both a 


936 


806 


24 


27 


5 


AovHarm 


Either b 


2118 


1719 


53 


375 


45 


AovHarm only c 


Either b 


450 


120 


4 


317 


79 


BLS 


EPD 


463 


79 


60 


39 


281 


BLS 


TFA 


444 


160 


81 


39 


160 


BLS 


Both a 


155 


62 


52 


4 


13 


BLS 


Either b 


752 


177 


89 


74 


428 


BLS only c 


Either b 


439 


l d 


21 


64 d 


428 


DACF 


EPD 


1491 


620 





534 


243 


DACF 


TFA 


1190 


507 





318 


95 


DACF 


Both a 


353 


203 





27 


21 


DACF 


Either b 


2328 


924 





825 


317 


DACF only c 


Either b 


1344 


373 





582 


528 


Totals^ 




4579 


2239 


95 


1176 


1105 



a Stars that arc selected by this method for both the EPD and TFA light curves. Note 
that stars that are automatically selected by this method for both EPD and TFA light 
curves but for which the by eye classification differs between the two light curve types will 
be included in the number of automatically selected variables, but will not be included in 
any of the subsequent columns for this row.^ Stars that arc selected by this method for 
cither the EPD or the TFA light curves. c Stars that arc selected exclusively by this method 
for either the EPD or the TFA light curves. For the Non-EB Var. and EB types stars are 
included in this row if they were only classified as that type during the visual inspection 
for this method. The questionable column lists the total number of stars that were flagged 
as questionable by this method and were either rejected during the visual inspection, or 
not selected, by all other methods. The rejected column lists the total number of stars that 
were rejected during the visual inspection for this method and did not pass the automatic 
selection for any of the other methods.** These stars are not included in the catalog of 
variables. For the EB and non-EB variable classes we list the total number of stars that 
were classified as this type during the visual inspection for at least one method. Note that 
36 stars arc flagged as EBs during the visual inspection for one method and as non-EB 
variables during the visual inspection for another method, so a total of 2298 stars are 
flagged as either an EB or as a robust non-EB variable for at least one method. This total 
does not include flare stars selected in § 16.31 unless they are also identified as a variable 
by AoV, AoVHarm, BLS or DACF. Including fiarc stars, the total number of stars with 
a robust variability detection is 2321. Of these, 1928 are not flagged as a probable blend 
or as having a problematic amplitude. For the questionable detections wc list the total 
number of stars that are classified as questionable for at least one method and are not 
classified as a robust detection for any of the methods. For the rejections wc list the total 
number of stars that were rejected during the visual inspection for all methods by which 
they were automatically selected. 



to determine the significance of the quad ratic fit 
relati ve to the constant function fit (see iLuptonl 
Il993f l. For each of the iV pca k sets we record 
the time-lag of the peak and its error from the 
quadratic fit (t p and crt p ), the peak DACF value 
and its error (y p and cry p ), and the false alarm prob- 
ability of the fit from the F-test (Pr p ). Note that 
this is not the false alarm probability of finding any 
connected set of points in the DACF that is well-fit 
by a quadratic in a random signal. In general, that 
false alarm probability will be higher. 

Starting from the peak with the shortest time lag 
t p , identify the first peak with Pr p < Pri; m i and 
Up > y\im or with Pr p < Pru m 2, choose this peak 
as the period for the star. If there is no such peak 
in the DACF, then the star is not selected as a 
variable by this method. We adopt Pru m i = 10 -4 ; 
PTn m 2 and j/iim are determined independently for 



each subsample from the simulations. For a given 
subsample we take Pry lm 2 and yy lm to be the fifth 
smallest and largest values respectively from the 
simulations (i.e. the 99.5 percentile values). 

In Figure [5] we plot the minimum Pr p found in each 
DACF against the maximum y p found and show the 
adopted cutoffs. We plot the results from the observed 
light curves and the simulations for four representative 
subsamples. Figure [7] shows an example of the DACF 
and selected peak for a periodic variable, and for one of 
the least significant detections that pass our selection. 
The un-phased light curves for each of these cases are 
also presented. 

We find significant peaks in the DACF for a total of 
1491 EPD light curves and 1190 TFA light curves. We in- 
spect these by eye to eliminate obvious false alarms (typ- 
ically cases where a light curve shows significant scatter 
on a few nights, often of several magnitudes or more), 
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Fig. 4. — Period vs. S/N from AoV for 4 representative light 
curve subsamples. On the left column we plot the observed values 
showing the light curves that pass the selection (dark filled points) 
and the light curves that do not pass the selection (grey filled 
points) separately, on the right column we plot the results from 
the Monte Carlo simulation for the corresponding subsample. The 
lines show the adopted S/N cut-off as a function of period. Note 
that in addition to the cut-off shown with the line, we also reject 
light curves for which the peak period is close to one sidereal day 
or a harmonic of one sidereal day. 

we also note whether or not the detection appears to 
be robust and whether or not the period determination 
is likely to be accurate. For the EPD light curves we 
consider 465 of the detections to be robust and with a 
correct period, 155 to be robust but at the wrong pe- 
riod (typically the detected period is an integer multi- 
ple of the likely true period, see for example figure [7]), 
534 to be non-robust, and we reject 243 of the detec- 
tions as clear false alarms. For the TFA light curves 
the respective numbers are 353, 154, 318, and 95 (see 
also table [3|) . The distinction between a robust and a 
non-robust detection is subjective, typically we consider 
a detection to be robust if the variability in the phased 
or unphased light curve is obvious to the eye and/or the 
DACF shows a set of clear regularly spaced peaks. There 
are some cases where the DACF shows clear regularly 
spaced peaks; however the scatter in the light curve ap- 
pears to correlate with the phase. We consider many of 
these cases to be non-robust (most are light curves from 
4K fields with periods near the lunar cycle). We include 
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Fig. 5. — Same as Figure [4] here we show the results for the 
AoVHarm period search. 

all non-rejected detections in the final catalog together 
with flags indicating the reliability of the detection and 
the reliability of the period. 

3.3. BLS - Search for Eclipses 

The BLS algorithm, primarily used in searches for 
transiting planets, detects periodic box-like dips in a light 
curve. This algorithm may be more sensitive to detached 
binaries with sharp-featured light curves than the other 
methods used. For our implementation of the BLS algo- 
rithm we search 10,000 frequency points within a period 
range of 0.1 to 20.0 days. To search for long period events 
where the eclipse duration may be only a very short frac- 
tion of the orbital period, we repeated the search at a 
higher frequency resolution using 100,000 points within 
a period range of 1.0 to 20.0 days. At each trial frequency 
we bin the phased light curve into 200 bins, and search 
over fractional eclipse durations ranging from 0.01 to 0.1 
in phase. Figure [8] shows the S/N vs the period for the 
EPD and TFA light curves. We select 752 stars with 
S/N > 10.0 and with a period not close to 1 sidereal day 
or a harmonic of a sidereal day as potentially eclipsing 
systems. We do not use the Monte Carlo simulations to 
set the selection thresholds for this method because dis- 
tinguishing between red noise and an eclipse signal by 
eye is less ambiguous than distinguishing between red 
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Fig. 6. — Minimum false alarm probability Pr p vs. maximum 
peak value y p found in each DACF for 4 representative light curve 
subsamples. On the left column we plot the observed values show- 
ing the light curves that pass the selection (dark filled points) and 
the light curves that do not pass the selection (grey filled points) 
separately, on the right column we plot the results from the Monte 
Carlo simulation for the corresponding subsample. The lines show 
the adopted Pr p /y p cut-off. Note that the plotted Pr p and y p 
values for a given DACF may not come from the same peak. In 
selecting peaks from the DACF we require that both the Pr p and 
the y p values for a given peak pass the selection. This plot is meant 
only to provide an approximate visualization of the selection. 

noise and general periodic or quasi-periodic variability. 
The selected light curves are inspected by eye to identify 
eclipsing systems. A total of 89 candidate EB systems 
are found in this manner, 8 are found in the EPD light 
curves only, 29 are found in the TEA light curves only, 
and 52 are found in both the EPD and TFA light curves 
(see also table [3]). 

3.4. Search for Flares 

As noted in the introduction, flaring is a common phe- 
nomenon among K and M dwarfs. While we were in- 
specting the light curves of candidate variable stars we 
noticed a number of stars showing significant flares. We 
therefore decided to conduct a systematic search for flare 
events in the light curves. Most optical stellar flares show 
a very steep rise typically lastin g from a few seconds to 
several minutes. iKrautten (|1996l ) notes that flares can be 
divided into two classes based on their decay times: "im- 
pulsive" flares have decay times of a few minutes, to a 



few tens of minutes, while "long-decay" flares have decay 
times of up to a few hours. Due to the 5-minute sampling 
of the HATNct light curves, flares of the former type will 
only affect one or two observations in a light curve, while 
flares of the latter type might affect tens of observations. 
In general it is very difficult to determine whether a given 
outlier in a light curve is due to a flare or bad photometry 
without inspecting the images from which an observation 
was obtained. This is impractical to do for tens of thou- 
sands of light curves when each light curve may contain 
tens to hundreds of outliers. While observations that are 
potentially corrupted are flagged, in practice the auto- 
mated routines that generate these flags do not catch all 
cases of bad photometry. We therefore do not attempt to 
identify individual "impulsive" flares in the light curves, 
and instead conduct a statistical study of the frequency 
of these flares (§ I6.3p . Long-decay flares, on the other 
hand, may be searched for in an automated fashion if a 
functional form for the decay is assumed (this is simi- 
lar to searching fo r microlensing events, see for example 
iNataf et al.l l2009f ) . To search for long-decay flares we 
used the following algorithm: 

1. Compute mo, the median magnitude of the light 
curve, and Sq the median deviation from the me- 
dian. 

2. Identify all sets of consecutive points with m — 
too < — 3<5o- Let to be the time of the brightest 
observation in a given set, and let N be the num- 
ber of consecutive points following and including to 
with to — mo < — 28a. We proceed with the set if 
N > 3. 



3. Use 



Leyenberg-Marquardt algorithm 
N points a 



the 

(jMarquardtJ |1963f) to fit to the 



function of the form: 
m(t) = -2.5 log 10 (Ae~ (t ~* n)/r + l) + mi (7) 

where A, r and TOi are the free parameters. Here A 
is the peak intensity of the flare relative to the non- 
flaring intensity, r is the decay timescale, and mi is 
the magnitude of the star before the flare. For the 
initial values we take mi = too, t = 0.02 days, and 
A = l0 -0 - 4 ( mp ~ 
at the peak. 



1, where m p is the magnitude 



4. Perform an F-test on the statistic given in eq. [6] 
where Xjv-3 m this case lS the \ 2 value from fitting 
eq. [7] If the false alarm probability is greater than 
1%, reject the candidate. If not, increase the num- 
ber of points by one and repeat step [31 Continue 
as long as the false alarm probability decreases. 

5. We reject any flare candidate for which there are 
at least two other candidate flares from light curves 
in the same field that occur within 0.1 days of the 
flare candidate. 

6. Let the number of points with t < to and to — t < 
0.05 days be N^ e { OTe and the number of points with 
t > ti and t — t\ < 0.05 days be iV a f ter . Here 
t\ is the time of the last observation included in 
the fit. Reject the candidate flare if A/before < 2 or 
iVafter < 2. Also reject the candidate flare if A < 0., 
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Fig. 7. — Top: Examples of the DACF for one of the highest significance detections of variability (left) and for one of the lowest 
significance detections (right). In each case the dark line shows the quadratic fit to the DACF used to determine the period of variation. 
Bottom: Un-phased light curves for the two stars. Note that the period for the star on the right is likely half the value determined by the 
peak identification algorithm. In this case the peak at P ~ 10 days did not pass the cut on Pr p , while the second peak at P ~ 20 days did. 
The P ~ 10 days signal is the top peak in the AoV and AoVHarm periodograms for this star, however the the S/N for both periodograms 
is below our selection threshold. 



A > 10.0, r < 0.001 days, r > 0.5 days, A < 
a a, t < a T , or if the false alarm probability from 
step S] is greater than 0.1%. Here a a and a T are 
the formal uncertainties on A and r respectively. 
The selection on A is used to reject numerous light 
curves with significant outliers which appear to be 
due to artifacts in the data rather than flares. 

We apply the above algorithm to the non-composite 
EPD light curves (i.e. light curves from each field are 
processed independently for stars with light curves from 
multiple fields) . The algorithm is applied both on the raw 
EPD light curves, and on EPD light curves that are high- 
pass filtered by subtracting from each point the median 
of all points that are within 0.1 day of that point. There 
are a total of 23, 589 stars with EPD light curves that 
are analyzed, we exclude from the analysis 4 stars from 
the full sample for which only cr-clipped TFA light curves 
are available. A total of 320 candidate flare events from 
281 stars are selected. These are inspected by eye to yield 
the final sample of 64 flare events from 60 stars. Figure [5] 
shows two examples of these large-amplitude, long-decay 
flares. The identified flares have peak intensities that 
range from A = 0.09 to A = 4.21 and decay time-scales 
that range from r = 4 minutes to r = 1.7 hours. 

3.5. Comparison of Selection Methods 

As seen in Table [3] more stars pass the automatic se- 
lections for AoV or for DACF than pass the automatic 
selections for AoVHarm. On the other hand, there are 



more robust detections found by AoVHarm than by the 
other methods. The latter result may be due, in part, 
to a bias toward sinusoidal signals in the by-eye selec- 
tion. However, taking it at face value, it appears that 
the AoVHarm method is a more robust period-finder for 
this sample of light curves than either the AoV or the 
DACF methods. For the 1075 TFA light curves that 
are classified as robust detections during the visual in- 
spections for both the AoV and AoVHarm methods, the 
S /N for the AoVHarm detection is greater than the S /N 
for the AoV detection in all but 12 cases. This is in 
line with the long-known fact that Fourier-based period- 
finding methods generally yield higher S/N detections 
than phase-binning m ethods, even i n cases where the sig- 
nal is non-sinusoidal (|Kovacslll980D . 

A direct comparison between the DACF and the 
AoVHarm or AoV methods is difficult since the selection 
for DACF is quite different from the selections for AoV or 
AoVHarm. However, it is apparent from table [3] that, at 
least for our selection thresholds, the DACF detections 
are generally found to be less reliable than the AoVHarm 
or AoV detections. Again there may be a bias in the 
by-eye selection against the types of variables that pass 
DACF, however it is also likely that the increased com- 
plexity in identifying a period from the auto-correlation 
function results in more false alarms than are generated 
by the periodogram-based methods. 

As expected, BLS appears to be the best technique for 
identifying eclipsing binaries; 89 out of the 95 potential 
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Fig. 8.— Period vs. S/N from BLS for EPD and TFA light 
curves. The dark points show light curves that pass the S/N > 10. 
cut and do not have a period near a sidereal day or one of its 
harmonics. The grey points show light curves that do not pass this 
cut. 



eludes flux contributions from many stars. Blending is a 
particularly significant issue in high stellar density fields 
near the Galactic plane. A star blended with a nearby 
variable star may be incorrectly identified as a variable 
based on its light curve. If the stars are separated by 
more than a pixel or two, it may be possible to distin- 
guish the real variable from the blend by comparing the 
amplitudes of their light curves. However, because pho- 
tometry is only obtained for stars down to a limiting 
magnitude (the value used varies from field to field), in 
many cases we do not have light curves for all the faint 
neighbors near a given candidate variable star, so we can- 
not easily determine which star is the true variable. In 
these cases we can still give an indication of whether or 
not a candidate is likely to be a blend by determining 
the expected flux contribution from all neighboring stars 
to the candidate's light curve. 

To determine whether or not a candidate is blended 
with a nearby variable star that has a light curve, we 
measure the peak-to-peak light curve amplitude (in flux) 
of all stars within 2' of the candidate. If any neighbor 
has an amplitude that is greater than twice the flux am- 
plitude of the candidate, the candidate is flagged as a 
probable blend. If any neighbor has an amplitude that 
is between half and twice the flux amplitude of the can- 
didate, the candidate is flagged as a potential blend. If 
any neighbor has an amplitude that is between 10% and 
half the flux amplitude of the candidate, the candidate 
is flagged as an unlikely blend. And finally we flag the 
candidate as a non-blend if all neighbors have amplitudes 
that are less than 10% that of the candidate. We deter- 
mine the amplitude of a light curve by fitting to it 10 
different Fourier series of the form: 



JV 



m(t) = mo + a* sin (2mt/P + <fii) 



(8) 
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Fig. 9. — Two examples of large-amplitude, long-decay flares 
seen in the HATNet light curves. In each case the solid line shows 
the fit of eq. [7] to the light curve that was done in automatically 
identifying candidate flares. 

eclipsing binaries identified in our survey are selected by 
BLS, including 21 that are identified exlusively by BLS. 
The second most successful method for identifying EBs 
is the AoV method which identified 70 of the 95 potential 
eclipsing binaries, and exclusively identified 4 of them. 

4. VARIABILITY BLENDING 

While the wide FOV of the HATNet telescopes allows 
a significant number of bright stars to be simultaneously 
observed, the downside to this design is that the pixel 
scale is necessarily large, so a given light curve often in- 



with N ranging from 1 to 10. Here P is the period of 
the light curve. We perform an F-test to determine the 
significance of each fit relative to fitting a constant func- 
tion to the light curve, and choose the amplitude of the 
Fourier series with the lowest false alarm probability. If 
the lowest false alarm probability is greater than 10% we 
set the amplitude to zero. We try all periods identified 
for each candidate by the variability searches described 
in § [31 and adopt the largest amplitude found. For can- 
didates that have light curves from multiple fields, or 
that have both ISM and AP reductions, we do the am- 
plitude test on each separate field/reduction and adopt 
the most significant blending flag found for the candi- 
date. If the amplitude of the candidate variable star is 
set to zero for a given field/reduction we do not use that 
field/reduction in determining the blending flag. If this 
is true for all fields/reductions we flag the candidate as 
problematic. We use the EPD light curves in doing this 
test. 

To determine whether or not a candidate is potentially 
blended with a nearby faint variable star that does not 
have a light curve, we compare the observed amplitude of 
the candidate to its expected amplitude if a neighboring 
star were variable with an intrinsic amplitude of 1.0 mag. 
We assume that an amplitude of 1.0 mag is roughly the 
maximum value that one might expect for a short period 
variable star. If the measured amplitude is less than the 
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expected amplitude then we flag the candidate as a po- 
tential blend, if it is greater than the expected amplitude 
and less than twice the expected amplitude we flag the 
candidate as an unlikely blend, and if it is greater than 
twice the expected amplitude, we flag it as a non-blend. 
The test is done for all stars within 2' of the candidate 
that do not have a light curve, and we adopt the most 
significant blending flag found for the candidate. To de- 
termine the expected amplitude of the candidate star 
induced by the neighbor, we note that a star with mag- 
nitude mi located near a variable star with magnitude 
m 2 and amplitude Am 2 > 0, has an expected light curve 
amplitude that is given by 



lPoimanskill2002f) 8 . and to the Super- WASP catalogue of 
eriodic variables c oincident with ROSAT X-ray sources 



Ami,AP = 2.51og 10 /ilO 



\-0.4mi 



+ / 2 10 



-0.4(m2-Am2) 



-2.5io gl0 [Aicr - 4 ™ 1 + / 2 io-°- 4m2 ] 

for the case of aperture photometry, and by 
AmijsM = 

2.51og 10 [/xKT - 4 '" 1 + / 2 l(T - 4m2 (10°' 4A " 12 
-2.51og 10 [Aicr - 4 ™ 1 ] 



(9) 



-1 



(10) 



for the case of image subtraction photometry. The two 
expressions differ because in the ISM pipeline photom- 
etry is done on difference images (only differential flux 
is summed in the aperture), whereas in the AP pipeline 
photometry is done directly on the science images (all 
stellar flux is summed in the aperture). Here /i :2 is 
the fraction of the flux from star 1(2) that falls within 
the aperture, and we use the catalog values (transformed 
from 2MASS) for mi and m 2 . To determine j\ and / 2 we 
integrate the intersection between the circular aperture 
and a Gaussian PSF which we assume to have a FWHM 
of 2 pixels (this is a typical effective "seeing" for both 
the 2K and 4K images), we do not consider pixelation 
effects in making this estimate. 

We compare the results from the two blending tests for 
each candidate, and adopt the most significant blending 
flag from among the two tests for the catalog. We do not 
run the test on the candidate flare stars, unless the star 
was selected as a variable by another method as well. Out 
of the 3474 stars that are in either the first or the second 
catalog, 936 are flagged as unblended, 451 are flagged 
as unlikely blends, 1397 are flagged as potential blends, 
399 are flagged as probable blends, and 291 are found to 
have problematic amplitudes (cases where the amplitude 
measuring algorithm failed for the star in question). 

5. MATCH TO OTHER CATALOGS 

5.1. Match to Other Variable Star Surveys 

We match all 3496 stars selected as potential vari- 
ables to the c ombined General C atalogue of Variable 
Stars (GCVS; IS^usLeialJ [20061) the New Cata logue 
of Suspected Variable Stars (NSV; lKholopovlll982fl and 
its su pplement (NSVS; iKazarovets. Samus' fc Durlevichl 
Il998h 7 . We also match to the ROTSE catalog of vari- 
able stars (jAkerlof et alJ 12000( 1 . to the ALL Sky Au- 
tomated Survey Catalogue of Variable Stars (ACVS; 

7 The GCVS, NSV and NSVS were o btained from 
I http : / / www . sai . msu.su/ groups /cluster / gcvs / gcvs / 1 on 2009 
April 7 



Norton et al.| [2007). In all cases we use a 2' matching ra- 
dius. We use a large matching radius to include matches 
to known variables that may be blended with stars in 
our sample. In total 77 of our candidate variables lie 
within 2' of a source in one of these catalogs, meaning 
that 3419 are new identifications. This includes 36 that 
match to a source in the GCVS, 4 that match to a source 
in the NSV, 7 that match to a source in the NSVS, 4 that 
match to a source in the ACVS, 8 that match to a source 
in the ROTSE catalog (4 of which are in the GCVS as 
well), and 23 that match to a Super- WASP source (2 of 
these are in their catalogue of previously identified vari- 
ables). Two of the 36 candidate variables that match to 
a source in the GCVS (HAT-215-0001451 and HAT-215- 
0001491) actually match to the same source, V1097 Tau, 
a weak emission-line T Tauri star. Both stars are flagged 
as probable blends in our catalog, in this case HAT-215- 
0001491 is the correct variable while HAT-215-0001451 
is the blend. 

We inspect each of the 77 candidates with a poten- 
tial match and find that the match is incorrect for 39 
of them and correct for 38. For 23 of the 39 incorrect 
matches the candidate variable is flagged as a probable 
blend in our catalog. In 8 cases the candidate variable 
is flagged as a potential blend, in 4 cases it is flagged as 
an unlikely blend, in 3 cases it is flagged as unblended, 
and in 1 case the amplitude is considered problematic. 
The match appears to be correct for four of the candi- 
date variables flagged as probable blends. In addition 
to HAT-215-0001491, the stars HAT-239-0000221 and 
HAT-239-0000513 both match correctly to sources in the 
GCVS. These stars form a common proper motion, low 
mass binary system. Both stars are flagged as proba- 
ble blends in our catalog. Each matches separately, and 
correctly, to a flare star in the GCVS (V0647 Her and 
V0639 Her respectively). For the NSV-matching prob- 
able blend candidate HAT-121-0003519 the match may 
be correct, though the positional uncertainty of the NSV 
source is high. A variable star classification is not avail- 
able for this source. 

The 16 variables that match correctly to a source in the 
GCVS include 4 BY Draconis-type rotational variables, 
6 UV Ceti-type flare stares, 1 INT class Orion variable 
of the T Tauri type, and 5 eclipsing systems. The EBs 
include the two W UMa-type contact systems DY CVn 
and V1104 Her, the two Algol-type systems DK CVn 
and V1001 Cas, and the M3V7 white dwarf EB DE CVn 
((van den Besselaar et alJl2007t ). 

Table [TH at the end of the paper, lists the first ten 
cross-identifications, the full table is available electroni- 
cally with the rest of the catalog. 

5.2. Match to ROSAT 

We match all 3496 stars selected as potent ial variables 
to th e ROSAT All-sky survey source catalog (jVoges et al.1 
1999) using the US National Virtual Observatory catalog 
matching utilities. We use a 3.5er positional matching 
criterion. A total of 248 of the variables match to an 
X-ray source, including 237 stars in our primary catalog 



8 Version IT obtained 

http://www. astrouw.edu. pl/asas/?page=catalogues 



from 
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of periodic variables, 14 of the EBs, and 24 of the 60 
flare stars. A few of the variable stars are close neighbors 
where one is likely to be a blend of the other, so there are 
243 distinct X-ray sources that are matched to. Table [T"7l 
at the end of the paper gives the cross-matches. The 
full table is available electronically with the rest of the 
catalog. In Section 16.21 we discuss the X-ray properties 
of the rotational variables. 

6. DISCUSSION 

6.1. Eclipsing Binaries 

The 95 stars that we identify as potential EBs have 
periods ranging from P = 0.193 days to P = 24.381 days. 
We flag 11 of the candidate EBs as probable blends, 23 as 
potential blends, 21 as unlikely blends, 35 as unblended, 
and 5 as having problematic amplitudes. Figure flOl shows 
phased light curves for 12 of the EBs. 
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Fig. 10. — Example phased light curves for 12 of the 95 potential 
EBs found in the survey. The period listed is in days. 



J04404+3127B, respectively, which form a com- 
mon proper-motion 15" binary system. The two 
stars are both selected as candidate EBs with the 
same period, and are both flagged as potential 
blends in the catalog; based on a visual inspection 
of the light curves we conclude that the fainter 
component (HAT-216-0003316) is most likely the 
true P = 2.048 day EB. The fainter component, 
which has spectral type M3 on SIMBAD, also 
matches to the X-ray source RX J0440. 3+3126. 

3. HAT-127-0008153: matches to CCDM 
J03041+4203B, which is the fainter compo- 
nent in a common proper-motion 20" binary 
system. This star is flagged as a potential blend, 
the brighter component appears to match to the 
X-ray source 1RXS J030403.8+420319. Based on 
a visual inspection of the light curves we conclude 
that HAT-127-0008153 is likely the true variable. 

4. HAT-169-0003847: is 24" from the Super-WASP 
variable 1SWASP J034433.95+395948.0, the two 
stars are blended in the HATNet images, however 
from a visual inspection of the light curves we con- 
clude that HAT-169-0003847 is likely the true vari- 
able. 

5. HAT- 192- 000 1841: is 46" from a high proper- 
motion, KO star BD+41 2679. The two stars may 
be members of a common proper-motion binary 
system (the former has a proper motion of 65.79, 
— 151.77 mas/yr in RA and DEC respectively, while 
the latter has 89.76, -117.14 mas/yr). 

6. HAT-169-0003847: is flagged as an unlikely blend, 
and we confirm that BD+41 2679 is not the true 
variable. 

7. HAT-193-0008020: matches to GSC 03063-02208 
and has a spectral type of M Oe listed on SIMBAD 
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8. HAT-216-0007033: matches to the X-ray source 
RX J0436.1+2733 and has spectral type M4 listed 
on SIMBAD. 

9. HAT-341-0019185: is 43" from TYC 1097-291-1, 
the two stars appear to be members of a common 
proper motion binary system. HAT-341-0019185 is 
flagged as a probable blend in our catalog, though 
it does not appear that TYC 1097-291-1 is the real 
variable. 



In addition to matching the candidate EBs to other 
variable star catalogs (§ 15. lj) and to ROSAT (§ 15. 2j) we 
also checked for matches to previously studied objects 
using SIMBAD. The following objects had noteworthy 
matches: 

1. HAT-148-0000574: matches to the X-ray source 
1RXS J154727.5+450803, and w as previously dis- 
covere d to be an SB2 system by iMochnacki et al.l 
(|2002f ). we discuss this system in detail in § 16.1.11 

2. HAT-216-0002918 and HAT-216-0003316: 
match to CCDM J04404+3127A and CCDM 



6.1.1. The Low-mass EB 1RXS J154727. 5+450803 

The EB HAT-148-0000574 matches to 1RXS 
J154727.5+450803. Using RV observations ob- 
tained with the Cassegrain spectrograph on the 
David Dunlap O b serva tory (DDO) 1.88 m telescope 9 , 
IMochnacki et al.l (|2002l ) found that this object is a 
P = 3.54997 ± 0.00005 day double-lined spectroscopic 
binary system with component masses > 0.26 Mq. 
This system, however, was not previously known to be 
eclipsing. Here we combine the published RV curves 

9 Based on data obtained at the David Dunlap Observatory, 
University of Toronto 
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from iMochnacki et all (|2002t ) with the HATNet I-band 
light curve to provide preliminary estimates for the 
masses and radii of the component stars. 

Figure [Til shows the EPD HATNet light curve phased 
at a period of P = 3.550018 days together with a model 
fit, while figure [12] sho ws a fit to the radial v elocity ob- 
servations taken from IMochnacki et al.l {2002). Note the 
out of eclipse variations in the light curve, presumably 
due to spots on one or both of the components, which 
indicates that the rotation period of one or both of the 
stars is tidally locked to the orbital period. Since the 
HATNet light curve is not of high enough quality to 
measure the radii to better than a few percent preci- 
sion, we do not fit a detailed spot model to the light 
curve, and instead simply fit a harmonic series to the 
out of eclipse observations and then subtract it from the 
full light curve. We model the light curve using the JK- 



TEBOP program ( 


Southworth. Maxted & Smallevll2004 


ISouthworth et al.1 


20041) which is based on the Eclips- 



ing Binary Orbit Program (EBOP ; iPopper fc Etzellll981l : 
lEtzellll981fc iNelson fe Davislll972f ). but includes more so- 
phisticated minimization a nd error analy sis routines. We 
used the DEBiL program (|Devor!l2005 (i to measure the 
eclipse minimum times from the light curve, which in 
turn were used with the RV curves to determine the 
ephemeris. In modelin g the RV curves we fix e = 
(jMochnacki et al.ll2002t found e = 0.008 ± 0.007), and 
we fix = R2/R1 — 1.0 in modelling the light curve 
given q — 1.00 ± 0.02 from the fit to the RV curves (the 
light curve is not precise enough to provide a meaningful 
constraint on fc). For completeness we note that we as- 
sumed quadratic limb da rkening coeffi cients of a — 0.257, 
b = 0.586 for both stars (|Claretll2000T ). which are appro- 
priate for a T e g = 3000 K, \ogg — 4.5, solar metallic- 
ity star. The results are insensitive to the adopted limb 
darkening coefficients; we also performed the fit using the 
coefficients appropriate for a T c g = 4000 K, log g = 4.5 
star and found negligible differences in the resulting pa- 
rameters and uncertainties. The parameters for the sys- 
tem are given in tableHJ Note that the la errors given on 
the masses and radii are determined from a Monte Carlo 
simulation (|Southworth et al]|2005h . These are likely to 
be overly optimistic given the inaccurate treatment of 
the spots, and our assumption that the component radii 
are equal. 

Table [5] lists the masses and radii of the 4 other known 
double-lined detached EBs with at least one main se- 
quence component that has a mass less than O.3M0. 
We do not include RR Cae which is a white-dwarf/M- 
dwarf EB that has presumably undergone mass transfer 
(|Maxtedll2007l ). In figure [TBI we plot the mass-radius re- 
lation for stars in the range 0.15M Q < M < 0.3M©. 
Like the components of CM Dra and the secondary of 
GJ 3236, and unlike the components of SDSS-MEB- 
1 and the secondary of 2MASSJ04463285+190432, the 
components of 1RXS J154727.5+4 50803 have radii that 
are larger than predicted from the iBaraffe et al.l (|1998l ) 
isochrones (if the age is > 200 Myr). The radii are ~ 10% 
larger than the predicted radius in the 1.0 Gyr, solar 
metallicity isochrone. High precision photometric and 
spectroscopic follow-up observations, and a more sophis- 
ticated analysis of the data are need to confirm this. 

6.2. Rotational Variables 



TABLE 4 

Parameters for the EB 1RXS J154727. 5+450803 



Coordinates and Photometry 


RA (J2000) 


15:47:27.42 a 


DEC (J2000) 


±45:07:51.39 a 


Proper Motions 


[mas/yr] -259, 200 b 


J 


9.082 mag c 


H 


8.463 mag c 


K 


8.215 mag c 


Ephemerides 


P 


3.5500184 ± 0.0000018 day 


HJD 


2451232.89534 ± 0.00094 


Physical Parameters 


Mi 


0.2576 ± 0.0085 M Q 


M 2 


0.2585 ± 0.0080 M Q 


Ri = R 2 


0.2895 ± 0.0068 R e 


RV Fit Parameters 


7 


-21.21 ±0.41 km/s 


Kt 


55.98 ± 0.76 km/s 


K 2 


55.78 ± 0.83 km/s 


e 


0.0 d 


LC Fit Parameters 


J2/J1 


1.0734 ±0.030 


(Ri + R 2 )/a 


0.0737 ± 0.0014 


R2/R1 


1.0 e 


i 


86.673° ± 0.068° 



a SIMBAD b PPMX C 2MASS d fixed fixed based on q = 
1.004 ± 0.020 from fitting the RV curves. 



Figure [JJ] shows example phased light curves for 12 
randomly selected rotational variables found in our sur- 
vey. For the following analysis we only consider stars in 
our variables catalog that are identified as reliable detec- 
tions for at least one search algorithm and that are not 
flagged as probable blends or as having problematic am- 
plitudes. In order of preference, we adopt the AoVHarm 
TFA, AoVHarm EPD, AoV TFA, AoV EPD, DACF TFA 
or DACF EPD period for the star, choosing the period 
found by the first method for which the detection is con- 
sidered secure. 

6.2.1. Period- Amplitude Relation 

For FGK stars there is a well-known anti-correlation 
between stellar activity measured from emission in the H 
and K line cores, or from Ha emission, and the Rossby 
number {Ro, the ratio of the rotation p eriod to the 
chara cteristic time scale of convection, see iNoves et al.l 
1984), which saturates for short periods. Similar anti- 
correlations with saturation have been seen between 
Rg and the X-ray to bolometric luminosity ratio (e.g. 
iPizzolato et al.l l2003f ) , and between Ro and the am- 
plitude of photometric va riability (e.g. iMessina et al.1 
120011 ; lHartman et al.l 120091 ). Main sequence stars with 
M < O.35M0 are fully convective (|Chabrier fc Baraffd 
1997), so one might expect that the rotation-activity 
relation breaks down, or significantly changes, at this 
mass. Despite this expectation, several studies have 
indicated that the rotation-activity relation (measured 
using v s m i and Ha respectively) continues fo r late 
M-dwarfs (iDelfosse et all 119981: iMohantv fc Basril 120031: 
iReiners fe Basrill2007D . Recently, however. IWest fe Basril 
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TABLE 5 

Other Double-Lined EBs with a very late M dwarf component 



Name 


Period [days] 


A*i [Mq] 


M 2 [Mq] 


i?l [Rq] 


R 2 [Rq] 


CM Dra 


1.27 


0.2310 ± 0.0009 


0.2141 ±0.0010 


0.2534 ±0.0019 


0.2396 ± 0.0015 


SDSS-MEB-1 


0.407 


0.272 ± 0.020 


0.240 ±0.022 


0.268 ±0.010 


0.248 ± 0.009 


2MASSJ04463285+190432 


0.619 


0.47 ±0.05 


0.19 ±0.02 


0.57 ±0.02 


0.21 ±0.01 


GJ 3236 


0.771 


0.376 ±0.016 


0.281 ±0.015 


0.3795 ± 0.0084 


0.300 ± 0.015 



References. — CM Dra: Morales et alj 120081) ; ILacvl 119771k IMetcalfe eFail (1199611 : SDSS-MEB-1: IBlake et aT] (120081) : 
2MASSJ04463285+190432: IHebb et al.l (120061k G J 3236: the parameters listed for this system are determined by giving 



equal weight to the three models in llrwin^eT al. (2009b) 




0.4 0.42 0.44 0.46 0.48 0.5 0.52 0.54 0.56 0.58 
Phase 

Fig. 11. — Phased HATNet I-band light curve for the low-mass 
EB 1RXS J154727.5+450803. The top panel shows the full EPD 
light curve, the bottom two panels show a model fit to the two 
eclipses after subtracting a harmonic series fit to the out of eclipse 
portion of the light curve. 



(2009) have found that rotation and activity may not al- 
ways be linked for these stars. 

In figure [T5] we plot the rotation period against the 
peak-to-peak amplitude for stars in several color bins. 
The peak-to-peak amplitude is calculated for the EPD 
light curves as described in section 01 for stars observed 
in multiple fields we take the amplitude of the combined 
light curve. Stars without an available EPD light curve, 
or for which the amplitude measuring algorithm failed, 
are not included in the plot. A total of 1525 stars are 
included in the plot. For stars with V — K$ < 5.0 (cor- 
responding roughly to M > 0.25M©) the photometric 
amplitude and the period are anti-correlated at high sig- 
nificance. There appears to be a cut-off period, such that 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
Phase 

Fig. 12. — Circular orbit fit to the RV curves for the low-mass 
EB 1RXS J154727.5+4 50803. The observations are taken from 
IMochnacki et al.l | |2U02T) . 



in 
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CO 
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0.16 0.18 



0.22 0.24 0.26 0.28 
Mass 



Fig. 13. — Mass-radius relation for 7 main sequence stars in 
double-lined DEBs with M < 0.3M Q . The two points with the 
smallest error bars are the components of CM Dra, the open cir- 
cles are the components of SDSS-MEB-1, the open square is the 
secondary component of 2MASSJ04463285±190432, the X is the 
secondary component of GJ 3236, and the filled square marks the 
two components of 1RXS J154727.5+450803. The solid lin e show s 
the 1.0 Gyr, solar metallicity isochronc from Baraffe ct al. (1998). 
Note that the error bars for 1RXS J154727.5+450803 do not incor- 
porate systematic errors that may result from not properly mod- 
elling the spots or allowing the stars to have unequal radii. 

the period and amplitude are uncorrelated for stars with 
periods shorter than the cut-off, and are anti-correlated 
for stars with periods longer than the cut-off. In other 
words, the rel ations are satur a ted b elow a critical ro- 
tation period. lHartman et alj (|2009f ) find a saturation 
threshold of Ro = 0.31, which for a 0.6 Mq star corre- 
sponds to a period of ~ 8 days (assuming (B — V)q = 1.32 
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Fig. 14. — Phased EPD light curves for 12 randomly selected 
rotational variables found in the survey. The period listed is in 
days. The zero-point phase is arbitrary. For this figure we use the 
period found by applying AoVHarm to the TFA light curves. 

for stars of this mass, and using the empirical relation 
between the conv ective time scale and {B — V)o from 
iNoves et al.l 11984) . This is consistent with what we 
find for the bluest stars in our sample. For stars with 

V - K s > 5.0 (M < 0.25 M ), the period and ampli- 
tude are not significantly correlated, at least for periods 
< 30 days. This result suggests that the distribution of 
starspots on late M dwarfs is uncorrelated with rotation 
period over a large period range, and is perhaps at odds 
with Ha/usini studies which indicate a drop in activity 
for very late M-dwarf s tars with vsini < 10 km/s (e.g. 
iMohantv k Basrll [20031 ). 

6.2.2. Period-Color/Period-Mass Relation 

In figure [H] we compare the distribution of V ~ Ks 
colors for periodic variables to the distribution for all 
stars in the sample. We also show the fraction of stars 
that are variable with peak-to-peak amplitudes greater 
than 0.01 mag as a function of V — Ks- 

The plotted relation has been corrected for complete- 
ness by conducting sinusoid injection/recovery simu- 
lations to estimate our detection efficiency. In con- 
ducting these tests we divide the sample into 90 pc- 
riod/amplitude/color bins. We use color bins of 2.0 < 

V - K s < 3.5, 3.5 < V - K s < 4.0, 4.0 < V - K s < 4.5, 
4.5 < V ~ K s < 5.0 and 5.0 < V - K s < 6.0, three 
period bins of 0.1-1 day, 1-10 days, and 10-100 days, 
and 5 amplitude bins logarithmically spanning 0.01 to 
1.0 mag. While ideally a much finer grid would be used 
for these simulations, we were limited by computational 
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Fig. 15. — Rotation period vs. peak-to-peak photometric ampli- 
tude for 1525 stars. We divide the sample into 4 bins by V — Ks 
color, corresponding roughly to M > 0.6Mq, 0.5Mq < M < 
0.6M Q , 0.25Mq < M < O.5M , and'M < 0.25M Q , from~blue to 
red. We also list the Spearman rank-order correlation coefficient 
and the statistical significance of the correlation for each sample 
(note that negative values of rg imply that the period and ampli- 
tude are anti-correlated) . For stars with V — Kg < 5.0 the period 
and peak-to-peak amplitude are anti-correlate with high signifi- 
cance. The relation appears to be saturated for periods < 5 days, 
with hints that the saturation period increases for decreasing stel- 
lar mass. For stars with V - K$ > 5.0 (M < O.25M ), the period 
and amplitude are not significantly correlated for P < 30 days. 



resources to this fairly coarse sampling. For each bin 
we randomly select 1000 stars with the appropriate color 
(for color bins with fewer than 1000 stars we select with 
replacement). For each selected star we then choose a 
random period and amplitude drawn from uniform-log 
distributions over the bin and inject a sine curve with 
that period/amplitude and a random phase into the light 
curve of the star. If both EPD and TFA light curves are 
available for the star we inject the same signal into both 
light curves. We do not reduce the amplitude for the 
injection into the TFA light curve, this may cause us 
to slightly overestimate our detection efficiency. If the 
star was identified as a variable or a potential variable 
by our survey, we first remove the true variable signal 
from the light curve by fitting a harmonic series to the 
phased light curve before injecting the simulated signal. 
We then process the simulated signals through the AoV, 
AoVHarm and DACF algorithms using the same selec- 
tion parameters as used for selecting the real variables. 
We do not apply by-eye selections on the simulated light 
curves, so our detection efficiency may be overestimated 
(particularly for longer period stars where the by-eye se- 
lection tended to be stricter than the automatic cuts). 
To get the completeness corrected variability fraction we 
weight each real detected variable by 1// where / is the 
fraction of simulated signals that are recovered for the pe- 
riod/amplitude/color bin that the real variable falls in. 
We find that we are roughly ~ 80% complete over our 
sample of stars for peak-to-peak amplitudes greater than 
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0.01 mag and periods between 0.1 and 100 days. Consid- 
ering the recovery fraction separately for each variability 
method, we find that ~ 91% of the simulations are recov- 
ered by AoVHarm, and in ~ 97% of the recovered sim- 
ulations the recovered frequency is within 0.001 day -1 
of the injected frequency. For AoV we again find that 
~ 91% of the simulations are recovered, but that frac- 
tion of recovered simulations where the recovered fre- 
quency is within 0.001 day~ x of the injected frequency is 
~ 93% in this case. For the DACF method we find that 
only ~ 34% of the simulations are recovered, although if 
we consider only simulations where the injected period 
is greater than 10 days the fraction of simulations that 
are recovered by DACF is then ~ 81%. The fraction of 
these latter recoveries for which the recovered frequency 
is within 0.001 day -1 of the injected frequency is only 
~ 51% in this case, however. For the simulations the 
recovery frequency is relatively insensitive to the period 
and color and depends most significantly on the ampli- 
tude. For amplitudes between 0.01 mag and 0.022 mag 
the recovery fraction is ~ 65%, whereas for amplitudes 
above 0.05 mag the recovery fraction is ~ 90%. Above 
0.05 mag the recovery fraction is independent of ampli- 
tude. As noted above, the estimated completeness of 
~ 80% is likely to be overly optimistic since we do not 
include the by-eye selection, do not account for the re- 
duction in signal amplitude by TFA and use relatively 
"easy-to-find" sinusoid signals. Nonetheless we expect 
that the recovery frequency is > 70%. 

As seen in figure [TH] the fraction of stars that are de- 
tected as variables increases steeply with decreasing stel- 
lar mass. While only ~ 3% of stars with M > 0.7 Mq 
are found to be variable, approximately half of the stars 
with M < 0.2 Mq are detected as variables with peak- 
to-peak amplitudes > 0.01 mag (fig. IPo]) . We find that 
an exponential relation of the form 

Var. Frac. = (0.0034 ± 0.0008) e (°- 84±0 06 '( y - Ks ' (11) 

fits the observed relation over the color range 2.0 < V — 
K s < 6.0. 

In figure [17] we show the relation between period and 
color. For stars with V — Ks < 4.5 the measured dis- 
tribution of log P is peaked at ~ 20 days with a broad 
tail toward shorter periods and a more rapid drop-off for 
longer periods. Note that the cut-off for longer periods 
may be due to the correlation between period and am- 
plitude for these stars; stars with periods longer than 
~ 20 days may be harder to detect and not intrinsically 
rare. The peak of the distribution appears to be corre- 
lated with color such that redder stars are found more 
commonly at longer periods than bluer stars. For stars 
with V — Ks > 4.5 the distribution changes significantly 
such that the logP distribution appears to be more or 
less flat between ~ 0.3 and ~ 10 days, while red stars 
with P > 10 days are uncommon. The morphology is 
broadly consistent with what has been seen from other 
surveys. 

To demonstrate this, in Figure [18] we compare the 
mass-period distribution for stars in our survey to the re- 
sults from other surveys of field stars and open clusters. 
We choose to use mass for the comparison rather than 
observed colors because a consistent set of colors is not 
available for all surveys. The masses for stars in our sur- 
vey are estimated from their V — Ks colors (see Fig. [T6|) . 
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Fig. 16.— Top: The distribution of V - K s colors for the 1849 
stars in our catalog of periodic variables that are flagged as be- 
ing secure detections for at least one detection method and are 
not flagged as probable blends or as having problematic ampli- 
tudes, compared to the distribution of V — Kg colors for all 27, 560 
stars in the sample. On the top axis we show the corresponding 
main sequence stellar masses determined by combining the empiri- 
cal V — Kg vs . Mjj main-sequence for stars in the Solar neighbor- 
hood given by Johnso n & Appa l|2009j) with the mass-Mx relation 
from the Baraffc ct al. (1998J) 4.5 Gyr, solar-me tallicity isochrone 
with L m i x = 1.0. We used the relations from Carpenter (2001) 
to convert the CIT K magnitudes from the isochrones into the 
2MASS system. The distribution for variable stars is biased to- 
ward redder V — Ks colors relative to the distribution for all stars. 
The decrease in the total number of stars in the sample red- ward of 

V — Kg ~ 3.5 is due to the V-band magnitude limit of the PPMX 
survey. Bottom: The completeness corrected fraction of stars that 
are variable with peak-to-peak R or Iq amplitude > 0.01 mag as a 
function of V — Ks ; this fraction increases exponentially with color 
(solid line, eq. [TTJ. While only ~ 3% of stars with M > 0.7 Mq 
are found to be variable at the > 1% level, approximately half of 
the stars with M < 0.2 Mq are variable at this level. 

We take data from the M ount Wilson and Vienna-KPNO 
(jStrassmeier et al.1l2000[) samples of field stars with rota- 
tion periods. F or the Mount Wilson sample we use the 
com pilation by | Barnesl (I2007T) . the original data come s 
from lBaliunas et al.l(|1996[) and from lNoves et alj (|1984h . 
For the Vienna-KPNO sample we only consider stars 
which are listed as luminosity class V. We estimate 
the masses for stars in these samples using the same 

V — Ks to mass conversion that we use for our own 
sample. The V and Ks magnitudes for these field stars 
are taken from SIMBAD where available. We also com- 
pare our sample to four open clusters with ages between 



K /M Dwarf Variability Survey 



19 



100 




V-K s [mag] 



1000 



100 




100 



Period [days] 



Fig. 17. — Top: Period vs. V — Kg for 1785 stars in our cata- 
log of periodic variables that are flagged as being secure detections 
for at least one detection method, are not flagged as having incor- 
rect period determinations for that method, and are not flagged as 
probable blends or as having problematic amplitudes. There ap- 
pears to be a paucity of stars with P > 10 days and V — Kg > 4.5 
or V — Kg < 3. Bottom: The distributions of periods are shown 
for four color bins. For the three bluest color bins there appears 
to be a correlation between period and color, such that the mode 
period is longer for redder stars. For stars with V — Kg > 4.5, 
on the other hand, the period distribution appears to be biased to 
shorter periods. Using a K-S test, we find that the probability that 
the stars in any two of the different color bins are drawn from the 
same distribution is less than 0.01% for all combinations except the 
for the combination of the two intermediate color bins. For that 
combination the probability is ~ 0.7%. 



100- 200 Myr, including M35 (~ 180 Mvn lMeibom et all 
120091 ). and three clusters observed by the MONITOR 
project: M50 (~ 130 Myr ; llrwin et alJ l2009ah. NGC 2516 
(~ 150 Myr: llrwin et all 12007ft. and M34 (~ 200 Myr; 
llrwin et alj 12006ft . Finally, we compare our sample to 
the two oldest open clusters with sign ificant samples 
of ro tation periods: M37 (~ 550 Myr: lHartman et al.l 
2009ft, and the Hyades (~ 625 Myr; iRadick et alj 119871 



19951: IProsser et alj[l995h . For the MONITOR clusters 



we use the mass estimates given in their papers, these 
are based on the Jc-m ass relation from the appropri- 
ate [Barafeleriai] (|1998ft isochrone for the age/metallicity 
of each cluster. For M35, M37 and the Hyades we use 
the mass estimates derived fr om the V, Jc-ma ss relations 
from the appropriate YREC (|An et alJl2007ft isochrones 
for each cluster. F o r M35 we only include 214 stars from 
the iMeibom et~aT1 (|2009ft catalog that lie near the main 



sequence in V, B and Ic, and we exclude any stars which 
have a proper motion membership probability less than 
80%, or an R V membership probab ility less than 80% as 
determined bv lMeibom et al.l (|2009f ). We expect that the 
stars in our sample, and in the other field star samples, 
have a range of ages, but on average will be older than 
the stars in the open clusters. 

The sample of stars with rotation periods presented 
here is substantially richer than is available for other sur- 
veys of field stars. This is especially the case for later 
spectral types. The Mt. Wilson and Vienna-KPNO sur- 
veys primarily targeted G and early K stars, so there is 
not much overlap in stellar mass between those samples 
and our sample. The few Mt. Wilson stars with esti- 
mated masses < 0.8 M do show an anti-correlation be- 
tween mass and period, and have periods that arc longer 
than the majority of stars in our sample. The Vienna- 
KPNO stars, on the other hand, have periods that clus- 
ter around ~ 10 days, which is closer to the mode of the 
period distribution for stars of comparable mass in our 
sample. Since the Vienna-KPNO stars were selected as 
showing spectroscopic evidence for chromospheric activ- 
ity before their periods were measured photometrically, 
this sample is presumably more biased to shorter pe- 
riod active stars than the Mt. Wilson sample. Given the 
correlation between photometric amplitude and rotation 
period, we would also expect our sample to be biased 
toward shorter period stars relative to the Mt. Wilson 
sample. When compared to the open cluster samples 
we see clear evidence for evolution in the rotation pe- 
riods of low-mass stars. Stars in the younger clusters 
have shorter periods at a given mass, on average, than 
stars in our sample. The discrepancy becomes more ap- 
parent for stars with M < 0.5 Mq for which the pe- 
riod and mass appear to be positively correlated in the 
young cluster samples while they are anti-correlated in 
our sample. Looking at the ~ 600 Myr clusters, again 
the periods of stars in our sample are longer at a given 
mass, on average, than the periods of the cluster stars, 
however in this case the mode of the period distribution 
for the cluster stars appears to be closer to the mode 
of the period distribution for our sample than it is for 
the younger clusters. The lowest mass stars in the older 
clusters also do not show as significant a correlation be- 
tween mass and period as do the lowest mass stars in 
the younger clusters. For stars with M < 0.3 Mq the 
available field star and older open cluster samples are 
too sparse to draw any conclusions from when compar- 
ing to our sample. For the younger clusters, we note 
that the distribution of periods for the lowest mass stars 
is even more strongly peaked toward short periods than 
it is in our sample. This suggests that these stars do 
lose angular momentum over time, despite not having a 
tachocline. A more detailed comparison of these data to 
models of stellar angular momentum evolution is beyond 
the scope of this paper. 

6.2.3. Period-X-ray Relation 

Figure [TH] shows the fraction of variables that match 
to an X-ray source as a function of period. This frac- 
tion is constant at ~ 22% for periods less than ~4 days, 
for longer periods the fraction that match es to an X-ray 
source decreases as ~ P~ 8 . Following lAgiieros et al.l 
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Comparison to open clusters with 100 Myr < t < 500 Myr 
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Fig. 18. — A comparison of the mass-period distribution for stars 
in our survey to the results from other surveys. See the text for a 
description of the data sources. For clarity we make the comparison 
separately for field stars, open clusters with 100 Myr < t < 500 Myr 
and for two open clusters with t ~ 600 Myr. The rotation periods of 
stars in our sample at a given mass appear to be longer, on average, 
than the periods of stars in the open clusters, this is true across all 
mass ranges covered by our survey. The rotation periods from the 
Vienna-KPNO survey appear to be comparable, at a given mass, 
to the periods of stars in our sample, while the periods from the 
Mt. Wilson survey appear to be generally longer than the periods 
from our survey. It is likely that our survey and the Vienna-KPNO 
surveys are biased toward high-activity, shorter period stars than 
the Mt. Wilson survey is. 

(2009) we calculate the ratio of X-ray to J-band flux via 
loginX/x/jj) = log 10 fx + 0.4 J + 6.30 (12) 



where 1 count s _1 in the 0.1-2.4 keV energy range is 
assumed to correspond on average to fx = 10~ n erg 
cm~ 2 s _1 . In figuredOlwe plot the flux ratio as a function 
of rotation period for samples of variables separated by 
their V — K$ color. The X-ray flux is anti-correlated 
with the rotation period for stars with M > O.25M , for 
stars with M < 0.25Af Q there is still a hint of an anti- 
correlation, though it is of low statistical significance (the 
false alarm probability is ~ 20%). This result is similar 
to what we found for the photometric amplitude-period 
relation. 
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Fig. 19. — The fraction of non-EB periodic variable stars that 
match to a ROSAT source as a function of rotation period. The 
errorbars for the 3 longest period bins show lcr upper-limits. For 
stars with rotation periods less than ~4 days the fraction that 
matches to an X-ray source is constant at ~ 22%, for longer periods 
the fraction decreases as ~ P~ o s . 



6.3. Flares 

In Section 13.41 we conducted a search for large- 
amplitude long-duration flares, finding only 64 events 
in 60 out of 23, 589 stars with EPD light curves that 
were analyzed. There are likely to be many more flare 
events that have been observed but which cannot be eas- 
ily distinguished from non-Gaussian noise in an auto- 
mated fashion. The presence of these flares may, how- 
ever, be identified statistically by looking for an excess of 
bright outliers relative to faint outliers in the light curves 
that correlates with other observables, such as the rota- 
tion period. For each light curve we determine the excess 
fraction of bright n-5 outliers via 



fn 



N n ,- - N n> . 
N tnt 



(13) 



where iV nj _ is the number of points in the light curve 
with to — too < —nS, N n ,+ is the number of points with 
to — too > nS, too is the median of the light curve, S is 
the median value of |to — too|, and there are iVtot points 
in the light curve. We do this for n = 3, 5 and 10. 
Before calculating /„ we high-pass filter the light curve 
by subtracting from each point the median of all points 
that are within 0.1 day of that point. 

It is a well known fact that the distribution of magni- 
tudes in the light curve of a faint star is skewed about 
the median toward faint values. As such we expect /„ to 
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Fig. 20.— The ratio of 0.1-2.4 keV X-ray flux to J-band near 
infrared flux vs. the rotation period for non-EB periodic variable 
stars that match to a ROSAT source. We divide the sample into the 
same 4 color bins used in fig. 1151 We also list the Spearman rank- 
order correlation coefficient and the statistical significance of the 
correlation for each sample. The X-ray flux is anti-correlated with 
the rotation period at high significance for stars with V — Kg < 5.0. 
For stars with V — Kg > 5.0 (M < O.25M ) there is still a hint of 
an anti-correlation, though it is of low statistical significance. 

be less than zero and to decrease for fainter stars. Fig- 
ure [21] shows the excess fraction of bright outliers as a 
function of the light curve scatter for the high-pass fil- 
tered light curves. We compare the observed relation to 
the relation obtained for a simulated set of light curves 
generated using Poisson noise for the flux from the star, 
the sky and the sky annulus. We simulate one light 
curve for each observed light curve, using the observed 
sky fluxes and extinctions to determine the expected sky 
flux and star flux for each point in each light curve. We 
find that the median relation between the excess fraction 
of bright outliers and the light curve scatter is consis- 
tent with the expected relation; however there is greater 
scatter about this relation than expected from our ideal- 
ized noise model. It is unclear whether other sources of 
noise, such as systematic variations due to blending, flat- 
fielding errors, pixelation effects, etc., will yield a skewed 
distribution of magnitude values, and if so, in which di- 
rection the magnitude distribution will be skewed. It 
is unlikely, however, that this would be correlated with 
parameters such as the rotation period. Therefore, if a 
correlation is observed between rotation period and the 
excess fraction of bright outliers, then it is likely to be 
physical. 

In figure[H]we plot the median excess fraction of bright 
5-<5 and 10-6 outliers as a function of period for stars in 
our variables catalog that are identified as reliable detec- 
tions for at least one search algorithm, and that are not 
flagged as probable blends or as having problematic am- 
plitudes. In order of preference, we adopt the AoVHarm 
TFA, AoVHarm EPD, AoV TFA, AoV EPD, DACF TFA 
or DACF EPD period for the star. There appears to be a 
slight anti-correlation between period and excess fraction 
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Fig. 21. — The excess fraction of bright 55 outliers vs. the light 
curve scatter <5, shown for the ~27,000 observed light curves (top) 
and for a simulated sample of light curves (bottom). The dark 
points with the error bars show the median values of /s for the 
observed light curves, the solid line shows a power-law fit to the 
median values for the observed light curve, while the dashed line 
shows a power-law fit to the median values for the simulations. 

of bright 5-6 outliers such that stars with short periods 
(< 10 days) have slightly more bright 5-r5 outliers than 
faint 5-6 outliers compared to the median, while stars 
with long periods (> 10 days) have slightly fewer com- 
pared to the median. A Spearman rank-order correlation 
test rejects the null hypothesis that the excess fraction 
is not anti-correlated with the period at 99.92% confi- 
dence. For 10-6 outliers, the excess fraction of bright 
outliers also appears to be anti-correlated with the pe- 
riod, though at lower significance (99.5%). 

In figure [HI we also compare the distribution of peri- 
ods for stars with detected large-amplitude long-duration 
flares to the distribution for stars without such flare de- 
tections. A total of 31 of the 60 flare stars have a robust 
period determination and are not flagged as a probable 
blend. The period detection frequency of ~ 50% for flare 
stars is significantly higher than that for all other stars 
(< 7%). This result is expected if stellar flaring is associ- 
ated with the large starspots that give rise to continuous 
photometric variations. As seen in figure |2"21 the distri- 
bution of periods for flare stars is concentrated toward 
shorter periods than the distribution for non-flarc stars. 
The longest period found for a flare star is 18.2 days 
whereas 31% of the non-flare stars with period determi- 
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nations have periods greater than 18.2 days. Conducting 
a K-S test we find that the probability that the two sam- 
ples are drawn from the same distribution is less than 
1CT 6 . 

7. CONCLUSION 

In this paper we have presented the results of a vari- 
ability survey conducted with HATNet of field K and M 
dwarfs selected by color and proper motion. We used 
a variety of variability selection techniques to identify 
periodic and quasi-periodic variables, and have also con- 
ducted a search for large amplitude, long-duration flare 
events. We conducted Monte Carlo simulations of light 
curves with realistic noise properties to aid in setting the 
selection thresholds. Out of a total sample of 27, 560 
stars we selected 3496 that show potential variability, in- 
cluding 95 that show eclipses in their light curves, and 60 
that show flares. We inspected all automatically selected 
light curves by eye, and flagged 2321 stars (including 
those with flares) as being secure variability detections. 
Because the HATNet images have poor spatial resolu- 
tion, variability blending is a significant problem. We 
therefore implemented an automated routine to classify 
selected non-flare variables as probable blends, potential 
blends, unlikely blends, unblended or as having problem- 
atic amplitudes. Altogether we found 1928 variables that 
are classified as secure detections and are not classified 
as probable blends or as having problematic amplitudes 
(cases where the best-fit Fourier series to the light curve 
has a flase alarm probability greater than 10%). This 
includes 79 stars that show eclipses in their light curves. 
We identified 64 flare events in 60 stars, 38 of these stars 
are also selected as potential periodic or quasi-periodic 
variables (37 are considered reliable detections, of which 
36 have reliable period determinations and 31 have reli- 
able period determinations and are not flagged as proba- 
ble blends are as having problematic amplitude determi- 
nations). We matched the sample of potential variables 
to other catalogs, and found that 77 lie within 2' of a pre- 
viously identified variable, while 3419 do not. Including 
only flare stars and variables that are classified as secure 
detections and are not classified as probable blends or as 
having problematic amplitudes, 43 (including 7 EBs) lie 
within 2' of a previously identified variable, so that 1885 
are new identifications. 

One of the eclipsing binaries that we identified is the 
previously known SB2 system 1RXS J154727.5+450803. 
By combining the published RV curves for the compo- 
nent stars with the HATNet / band light curve, we ob- 
tained initial estimates for the masses and radii of the 
component stars (Tab. H|). The system is one of only 
a handful of known double-lined eclipsing binaries with 
component masses less than 0.3 Mq. While we caution 
that the errors on the component radii are likely to be 
underestimated due to systematic errors that have not 
been considered in this preliminary analysis, it is inter- 
esting that the radii do appear to be larger than predicted 
if the system is older than ~ 200 Myr. With a magni- 
tude V ~ 13.4, this system is only slightly fainter than 
the well-studied binary CM Dra (V ~ 12.90) which has 
been the anchor of the empirical mass-radius relation for 
very late M dwarfs. 1RXS J154727.5+450803 is thus a 
promising target for more detailed follow-up to obtain 
high precision measurements of the fundamental param- 



0.0025 



0.002 



0.0015 



0.001 



0.0005 



-0.0005 



-0.001 



-0.0015 



-0.002 



5-8 Outliers 
10-5 Outliers 




0.1 



1 10 
Period [days] 



100 



1400 



700 



5-8 Outliers 
10-8 Outliers 




200 400 600 800 1000 1200 1400 1600 1800 2000 
Rank Period 



1000 



100 




100 



Period [days] 



Fig. 22. — Top: The median excess fraction of bright 5(5 out- 
liers and 105 outliers vs. period for stars in our variables catalog 
that are identified as reliable detections for at least one search al- 
gorithm, and that are not flagged as probable blends or as having 
problematic amplitudes. In computing the excess fraction for each 
light curve we su btra ct the median excess fraction for the light 
curve scatter (fig. |2H . Center: The median rank of the excess 
fraction of bright 55 and 105 outliers vs. the rank period. The 
anti-correlation between period and the excess fraction of bright 
outliers is more apparent when the ranks, rather than the values, 
are plotted against each other. Bottom: Comparison between the 
period distributions for stars with detected large-amplitude long- 
duration flares, and stars without such a flare detected. Note the 
absence of long-period stars with flares. For clarity we have multi- 
plied the flare star period distribution by a factor of 10. 

eters of the component stars. With additional follow-up, 
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Fig. 23. — Top: The distribution of V— Kg colors for the 60 stars 
with high amplitude flares detected compared to the distribution of 
V — Kg colors for all 27, 560 stars in the sample. The distribution 
for flares stars is biased toward redder V — Kg colors relative to the 
distribution for all stars. Note that we use a 5 times higher binning 
resolution for the full sample. Bottom: The fraction of stars with a 
high amplitude flare detection is plotted against V — Kg . For stars 
with V — Kg < 4.0 less than ~ 0.1% of stars had a flare detected, 
for stars with V — Kg > 4.5 the fraction is > 1%. 



the large sample of 79 probable late-type eclipsing bi- 
naries presented in this paper should prove fruitful for 
further investigations of the fundamental parameters of 
low-mass stars. 

The majority of the variable stars that we have iden- 
tified are likely to be BY Dra type variables, with the 
measured period corresponding to the rotation period of 
the star. This is the largest sample of rotation periods 
presented to date for late-type field stars. We discussed 
a number of broad trends seen in the data, including an 
anti-correlation between the rotation period and the pho- 
tometric amplitude of variability, an exponential relation 
between V — Ks color and the fraction of stars that are 
variable, a positive correlation between period and the 
V — Ks color for stars with V — Ks < 4.5, a relative 
absence of stars with P > 10.0 days and V — Ks > 4.5, 
and an anti-correlation between the rotation period and 
the ratio of X-ray to J band flux. The correlations be- 
tween period and activity indicators including the ampli- 
tude of photometric variability and the X-ray emission 
are consistent with the well-known rotation-age-activity- 
mass relations for F, G, K and early M dwarfs. The 



data presented here may help in further refining these 
relations. Our data hints at a change in the rotation- 
activity connection for the least massive stars in the sam- 
ple (M < 0.25 Mq). The anti-correlation between period 
and amplitude appears to break down for these stars, and 
similarly the period- X-ray anti-correlation is less signifi- 
cant for these stars than for more massive stars. This is 
potentially at odds with previous studies which used Ha 
to trace activity and v sin i to infer rotation period, and 
found that the period-activity anti-correlation extends to 
very late- type M dwarfs. Comparing our sample to other 
field and open cluster samples, we find that the rotation 
periods of stars in our sample are generally longer than 
the periods found in open clusters with t < 620 Myr, 
which implies that K and M dwarf stars continue to lose 
angular momentum past the age of the Hyades. This 
appears to be true as well for stars with M < 0.25 M , 
though these stars generally have shorter periods than 
more massive stars in our sample. 

Finally we have conducted a search for flare events in 
our light curves, identifying 64 events in 60 stars. Due 
to the difficulty between distinguishing a flare from bad 
photometry in an automated way, there are likely to be 
many flare events in the light curves that we do not 
identify. We therefore do not attempt to draw conclu- 
sions about the total occurrence rate of flaring (for a re- 
cent determination of this fre quency using data fro m the 
Sloan Digital Sky Survey see lKowalski et all 120091 ) . We 
find that the distribution of V — Ks colors for flare stars 
is biased toward red colors, implying that the flare fre- 
quency increases with decrea sing stellar mass, which has 
been kn own for a long time ()Ambartsumvan et alJll970L 
see also lKowalski et aTll2009l ). We find that roughly half 
the flare stars are detected as periodic variables, which 
is a significantly higher fraction than for the full sam- 
ple of stars. This is in line with the expectation that 
stellar flaring is associated with the presence of signif- 
icant star spots , and is consistent with the finding by 
iKowalski et al.l (|2009t ) that the flaring frequency of ac- 
tive M dwarfs showing Ha emission is ~ 30 times higher 
than the flaring frequency of inactive M dwarfs. We also 
find that the distribution of periods for flare stars is 
biased toward shorter periods, again as expected from 
the rotation-activity connection. Finally we attempt to 
statistically identify flares by searching for excess bright 
outliers relative to faint outliers in the light curves. This 
excess appears to be anti-correlated with rotation period 
and provides further evidence that flares are more com- 
mon among rapidly rotating K and M dwarfs than among 
slower rotators. 
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APPENDIX 

A. MONTE CARLO SIMULATIONS TO DETERMINE SELECTION CRITERIA 

Here we describe the Monte Carlo simulations that we conducted to inform our choice of selection criteria for several 
of the variability detection algorithms that we discussed in § [3] The goal of these simulations is to establish thresholds 
on each variability statistic such that one may be fairly confident that light curves with values above the threshold 
exhibit real variability. Under the ideal case where the light curves consist of uniformly sampled Gaussian white noise, 
the probability distribution for many of the variability statistics may be calculated analytically, and a formal false 
alarm probability may be determined. For nonuniform time sampling, however, the analytic expressions generally 
do not hold. In this case one may still use an analytic approximation for the false alarm probability by including a 
parameter for the effective number of ind ependent trial frequenci es. This parameter must either be determined by 
conducting Monte Carlo simulations fe.g. iHorne fc Ba liunas 1986), or it can be estimated from the total number of 
data points and the total number of frequencies sampled (see lPress et al.l ll992). In general, however, the noise in a 
light curve is neither Gaussian nor uncorrelated in time. Temporally-correlated, low- fr equency ("red") noise appears 
to be ubiquitous in ground based variability surveys (e.g. iPont. Zucker fc Queloal2006l ): additionally the distribution 
of magnitude values in ground based light curves often show broad non-Gaussian wings. As we will show below, the 
HATNet light curves exhibit both of these effects. In both cases the analytic expressions derived for Gaussian white 
noise will underestimate the actual false alarm probability. In this case one must conduct Monte Carlo simulations of 
light curves with realistic noise properties to establish the false alarm probabilities. 

The formally proper way to proceed would be to conduct many thousand simulations for each real light curve to 
empirically determine the distribution of AoV, AoVHarm, etc. values at each frequency in the periodogram of the light 
curve. These distributions would then be used in selecting the peak frequency (greatest outlier) in the periodogram, 
and to determine the false alarm probability of finding that peak. Such a procedure is unpractical, however, for 
more than a few tens of stars given the non-negligible amount of time required to calculate periodograms for a single 
simulated light curve (calculating the AoVHarm periodogram for a typical HATNet light curve takes ~ 2 minutes 
running as a single process on a 2.5 GHz AMD 9850 Quad-Core Processor). Instead we opt for a simpler procedure of 
generating an ensemble of light curves meant to represent the observed sample of light curves if they consisted only 
of noise. We run the variability selection routines on the ensemble and compare the resulting period vs. S/N etc. 
distributions to the real distributions. Selection criteria are adopted such that few of the light curves in the simulated 
sample would pass the cuts. Because the noise properties of the light curves do depend on several parameters, we 
divide the full sample of light curves into several sub-samples and establish separate criteria for each subsample as 
described below. 

In general the noise properties of a light curve will depend on many factors including the brightness of the star, its 
proximity to other stars, and its location on the image, to name just a few. Exploring all of these factors, however, 
is not practical given the computational constraints. Here we make the assumption that the noise depends only on 
the light curve type (EPD or TFA), the chip format (2K or 4K CCDs), and the magnitude of the star. For each chip 
format and filtering method we divide the sample of light curves into magnitude bins. We use a binsize of 1.0 mag 
and ranges from 6.0 < I < 14.0 for the 2K CCDs and 7.0 < R < 15.0 for the 4K CCDs. For the moment we exclude 
composite light curves that include both 2K and 4K observations. We conduct the simulations separately for each bin. 

We discuss first how we generate the simulated light curves, we then discuss how the noise properties are determined 
for each sub-sample of light curves. To simulate a light curve with red noise we use the following procedure: 

1. Choose an observed light curve as a basis for the time sampling of the simulated light curve (the simulated light 
curve will be interpolated to this time basis at the last step). 

2. Generate a series of uncorrelated random magnitude measurements drawn from an empirical distribution of 
magnitudes p(m). We discuss below how this empirical distribution is determined. We will refer to this series 
as s(t). The total number of generated points N s should be a power of 2 (i.e. log 2 N s is an integer) and such 
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that NsA. > T where A is the minimum time separation between points in the observed light curve and T is the 
total time-span covered by the observed light curve. 

3. Use the fast Fourier transform (FFT) algorithm (e.g. lCoolev fc Tu kcv 1965) to calculate s{v), the discrete Fourier 
transform (DFT) of s(i). 



4. To add red noise to the simulated signal we apply the filter H(v) to s(v) to get r(v) 

f{v) = H{v)s(v). (Al) 

The filter H{v) satisfies 

H(v) = yfPtv) (A2) 

where P{v) is the power spectrum of the red noise. We discuss below how it is determined. Note that we assume 
H{v) is real valued. 



5. Apply the inverse FFT to f(y) to obtain r(t), a uniformly sampled time-series with the desired noise properties. 

6. Interpolate r(t) at the time values of the observed light curve to yield the simulated light curve in its final form. 
Note that because the variability selection algorithms are insensitive to scale changes in the magnitudes, we do 
not need to adjust the RMS of the simulated light curve to equal that of an observed light curve. 



To determine the distribution of magnitude values we divide each light curve by its standard deviation calculated 
with 5cr iterative clipping about the median. The normalized magnitude values from all the light curves in a sub- 
sample are then grouped into a single list, which we take as p(rn) in step above. Note that we do not attempt to 
remove variable stars, or to whiten the light curves before generating this empirical distribution. Figure [24] shows the 
histogram of normalized magnitude values for 2K, EPD light curves with 9.0 < I < 10.0. Note that the wings of the 
distribution are non-Gaussian. We find that a function of the form 

p(m) = ^exp(-m 2 /2cr, 2 „) + B cxp(- H°- 5 /C) (A3) 

fits the observed distribution, with a m = 0.886, C = 0.427 and A/B = 1.21 for the aforementioned sub-sample. In 
simulating the light curves we draw points from the empirical distribution p{m), we give the functional form of the 
distribution for reference only. As a test on our assumption that the wings of the distribution are not dominated by 
contributions from variable stars or stars with significant trends, we also computed the distribution for the sub-sample 
shown in figure [24] after excluding stars identified as potential variables in § [3] and median filtering the light curves 
(subtracting from each point in the light curve the median value of all points within 0.1 day of the given point); the 
differences between the distributions are negligible. 

To determine the power spectrum of the red noise, P{v), used in eq, IA2l we calculate the DFT power spectra for 100 
randomly selected light curves in each sub-sample. The power spectra are interpolated to a common frequency grid 
and scaled to have unit median. These power spectra are then median combined to yield the median power spectrum 
-Pmed(^)- Figure [231 shows examples of P m ed(v) for several of the light curve sub-samples. Note the general rise in the 
power spectra at lower frequencies as well as the excess power near the sidereal frequency ^ s ;d = 1.0027379 day - and 
its harmonics. At exactly z/ s i<j the power spectra drop to low values due to the decorrelations in the EPD routine. Also 
note that the power spectra of fainter stars, and of the TFA light curves show less frequency dependence (i.e. they are 
closer to a horizontal line, which is what one would see for white noise). As shown in figure [231 we fit a function of 
the form 

to P me d( I/ )j where a, b, c, d, e and / are fitting parameters which are listed for each sub-sample in table [SI The 
parameters a and e tend to decrease with increasing magnitude, and are generally smaller for the TFA sub-samples 
than for the EPD sub-samples. Otherwise, there are no clear trends in the parameter values. We use this analytic 
expression for P{v) in simulating the light curves (eq. IA2[) . 

Figure [26] shows examples of several simulated light curves together with examples of some real light curves that are 
not selected as variable. The general noise properties of the observed light curves, including an abundance of outliers, 
and the presence of red noise, are reproduced in the simulated light curves. 

It is important to note that there are some properties of the noise in the real light curves that are not modeled in 
the simulated light curves. In the real light curves the scatter of the magnitude measurements about the expected 
value varies in time. On nights with poor weather conditions the scatter is often higher than on photometric nights. 
In our simulated light curves the scatter is constant though the expected value may vary in time. Also note that if the 
red noise properties of the light curves vary significantly within a given sample, it is likely the median power spectrum 
will show less red noise than any given light curve. For this reason we adopt conservative selection criteria that pass 
well above the simulated variability distributions in the selections described in § [3] 
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Fig. 24. — The histogram of normalized A magnitude values for 2K, EPD light curves with 9.0 < I < 10.0 (light filled circles). The 
light curves that contribute to this distribution have been divided by th eir s tandard deviations calculated with 5cr iterative clipping about 
the median. The solid line shows a fit to the observed distribution (eg I A3t . Note that we use a logarithmic scale for the vertical axis to 
highlight the wings of the distribution which are non-Gaussian. 
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Fig. 25. — Median DFT power spectra for 4 representative sub-samples of light curves. In all cases there is a rise in power toward 
lower frequencies which is symptomatic of red noise. Light curves also generally show increased power near the sidereal frequency and its 
harmonics. At exactly the sidereal frequency the power drops due to the decorrelation against variations at this frequency when applying 
EPD. The red noise is less significant for fainter stars than for brighter stars, and it is less significant after trend filtering with TFA (note 
that this is general TFA procedure, not the signal- reconstruction mode). That some red noise remains in the post-TFA light curves suggests 
that for many of the fields the TFA pr o cedure was run with fewer than the ideal number of templates required to achieve optimum noise 
filt ering (Kovacs, Bakos & Noyes 2005; Szulagyi, Kovacs & Welch 2009). The solid lines show fits to the observed power spectra using 
eq. IA4I These fits are used to simulate light curves with red noise. Note that power is plotted on a logarithmic scale. 
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TABLE 6 

Parameters from fitting eo. IA4I to the median power spectra of 
different light curve sub-samples. 



Chip 


LC Type 


Mag. 


a 


b 


c 


d 


c 


f 


2K 


EPD 


6.5 


12.466 


1.897 


2.190 


0.179 


2.338 


0.281 


2K 


EPD 


7.5 


21.315 


1.954 


2.464 


0.487 


5.235 


0.273 


2K 


EPD 


8.5 


27.054 


1.910 


2.194 


0.875 


8.467 


0.195 


2K 


EPD 


9.5 


13.549 


2.054 


2.319 


1.114 


4.635 


0.136 


2K 


EPD 


10.5 


8.118 


2.184 


1.957 


0.265 


4.009 


0.437 


2K 


EPD 


11.5 


6.415 


2.488 


2.043 


0.190 


4.177 


0.508 


2K 


EPD 


12.5 


5.962 


2.540 


2.270 


0.258 


4.630 


0.584 


2K 


EPD 


13.5 


2.862 


2.842 


2.073 


0.191 


1.882 


0.557 


2K 


TEA 


6.5 


2.327 


2.899 


2.375 


0.600 


1.343 


0.336 


2K 


TEA 


7.5 


2.704 


2.875 


2.091 


0.713 


1.252 


0.184 


2K 


TEA 


8.5 


4.242 


2.683 


2.163 


1.001 


1.758 


0.030 


2K 


TEA 


9.5 


1.530 


2.793 


2.181 


0.962 


0.539 


0.000 


2K 


TEA 


10.5 


1.095 


2.794 


2.173 


0.933 


0.305 


0.003 


2K 


TEA 


11.5 


0.645 


3.315 


2.357 


0.899 


0.107 


0.000 


2K 


TEA 


12.5 


0.372 


3.763 


2.459 


0.949 


0.050 


0.000 


2K 


TEA 


13.5 


0.296 


3.828 


2.472 


0.958 


0.027 


0.000 


IK 


EPD 


7.5 


7.088 


2.454 


2.168 


0.001 


1.075 


0.281 


IK 


EPD 


8.5 


8.024 


2.483 


2.438 


0.006 


1.733 


0.385 


4K 


EPD 


9.5 


17.608 


2.211 


2.735 


0.915 


6.401 


0.297 


4K 


EPD 


10.5 


21.052 


2.081 


2.515 


0.767 


6.782 


0.254 


IK 


EPD 


11.5 


7.768 


2.361 


2.519 


0.415 


5.382 


0.570 


4K 


EPD 


12.5 


6.790 


2.257 


2.191 


0.167 


5.061 


0.641 


4K 


EPD 


13.5 


3.331 


2.748 


2.358 


0.246 


3.157 


0.663 


IK 


EPD 


14.5 


1.794 


3.330 


2.239 


0.121 


2.188 


0.738 


4K 


TEA 


7.5 


1.723 


3.507 


2.514 


0.022 


1.836 


0.540 


IK 


TEA 


8.5 


1.848 


3.459 


2.358 


0.767 


0.574 


0.064 


IK 


TEA 


9.5 


3.031 


2.868 


2.211 


0.676 


0.949 


0.132 


IK 


TEA 


10.5 


2.226 


2.670 


2.158 


0.791 


0.419 


0.000 


4K 


TEA 


11.5 


0.723 


3.361 


2.229 


0.909 


0.125 


0.000 


4K 


TEA 


12.5 


0.307 


3.857 


2.694 


0.998 


0.047 


0.000 


IK 


TEA 


13.5 


0.182 


3.908 


2.735 


1.041 


0.040 


0.001 


IK 


TEA 


14.5 


0.118 


4.340 


2.896 


1.084 


0.069 


0.000 



B. CATALOG OF VARIABLE STARS 

We combine the results from the various selection algorithms to form 3 catalogs of variable stars: those selected 
by the AoV, AoVHarm or DACF methods applied to either the EPD or TFA light curves and not flagged as an EB 
during the visual inspection for that method, a catalog of stars with eclipses in their light curves, and a catalog of 
flares. All three catalogs will be made available in electronic form from the CDS archive server 10 as well as from the 
HATNet website 11 . Light curves for the variable stars will be made available through the Harvard University Time 
Series Center 12 and through NStED 13 . There are a total of 3418 stars that are in the first catalog, 95 candidate EB 
stars in the second, and 64 flare events from 60 stars in the third. Note that some of the EBs and flare stars are 
included in the primary catalog as well. Altogether 3496 stars are included in at least one of the catalogs. 

In the primary catalog of periodic variables we include the internal HAT-ID of each source, its J2000 position from 
2MASS, the periods, S/N, peak-to-peak amplitudes and quality flags for the AoV and AoVHarm algorithms applied 
to the TFA and EPD light curves, the periods, period uncertainties, DACF peak heights, false alarm probabilities, 
peak-to-peak amplitude and quality flags for the DACF algorithm applied to the TFA and EPD light curves, and the 
blending flag described in §0J For convenience we also include the PPMX proper motion, the 2MASS photometry, and 
the estimated V magnitude for each source. For each selection method we list the peak-to-peak amplitude determined 
for both the EPD and TFA light curves phased at that period. See § 0] for a description of how the amplitude is 
measured. Note that the TFA light curves used for the amplitude determination have not had been processed with 
signal-reconstruction mode TFA, as a result these amplitudes are generally lower than the real amplitude of the star, 
and are lower than the amplitudes measured on the EPD light curves. Tables [71ITT1 at the end of the paper show the 
first ten entries in this catalog. For display purposes we split the catalog into 5 separate tables, in the full electronic 
version of the catalog the tables are joined. 

In the catalog of candidate eclipsing binary stars we include the HAT-ID, position, proper motion, 2MASS photom- 
etry, and V magnitude, together with the period, period uncertainty, epochs of primary and secondary eclipse (if the 
secondary is detected) and their respective uncertainties, the durations of the primary and secondary eclipses and the 
blending flag. Tables [T2"l[T3l at the end of the paper show the first ten entries in this catalog. 

In the catalog of flare events we include the HAT-ID, position, proper motion, 2MASS photometry, and V magnitude, 

10 http://cdsarc.u-strasbg.fr 

11 http://hatnet.hu/ 

12 http:/ /timemachine. iic.harvard.edu/ 

13 http://nsted.ipac.caltech.edu 
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Fig. 26. — Examples of simulated light curves (right column) compared with the observed light curves used as bases for the time sampling 
(left column). We show an example with significant red noise (top row), and one with very little red noise (bottom row). The observed 
light curves have been normalized to have unit standard deviation after an iterative 5cr clipping. 

together with the time of the flare peak, the peak intensity relative to the non- flaring intensity (A in Eq. [7]), and the 
decay time of the flare (r in Eq. [7]). Tables [TUfTSl at the end of the paper show the first ten entries in this catalog. 
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TABLE 7 

Catalog of candidate periodic variable stars. I: Astrometric and Photometric Data 





ID 


RA a 


DE a 


pmRA b 


pmDE b 


Jmag c 


Hmag c 


Kmag c 


Vmag 


Vref 1 


Bl c 








h:m:s 


deg:m:s 


mas/yr 


mas/yr 


mag 


mag 


mag 


mag 






HAT- 


■085- 


-0004278 


23:22:00.71 


+56:59:19.9 


392.08 


74.00 


9.473 


8.868 


8.655 


13.237 





2 


HAT- 


-085- 


-0053487 


23:21:24.35 


+56:25:14.3 


110.04 


-73.81 


12.350 


11.679 


11.566 


14.636 


2 


2 


HAT- 


-086- 


-0001701 


00:13:38.42 


+52:45:05.0 


35.13 


-36.72 


9.519 


8.763 


8.658 


12.136 


1 


1 


HAT- 


-086- 


■0002357 


23:51:38.51 


+49:45:21.3 


324.91 


290.31 


9.062 


8.400 


8.210 


11.437 





1 


HAT- 


-086- 


-0003825 


23:54:10.49 


+51:41:09.9 


249.11 


51.86 


9.574 


9.000 


8.745 


13.139 


1 


2 


HAT- 


■086- 


■0004017 


23:54:23.85 


+54:41:26.9 


35.30 


35.51 


9.565 


8.804 


8.643 


12.235 


1 


2 


HAT- 


-086- 


-0004248 


00:07:19.34 


+52:41:31.6 


325.06 


-3.32 


9.736 


9.157 


8.940 


12.568 





2 


HAT- 


■086- 


-0005153 


23:51:18.20 


+54:35:16.2 


76.13 


-7.77 


9.906 


9.239 


9.060 


12.822 


1 


2 


HAT- 


■086- 


-0005542 


00:16:56.01 


+50:37:49.0 


146.68 


-27.87 


9.980 


9.349 


9.137 


12.767 


1 


2 


HAT- 


-086- 


■0005602 


00:16:01.61 


+48:55:56.9 


-83.81 


51.47 


10.008 


9.409 


9.186 


12.870 


1 


3 



a J2000, from 2MASS b from PPMX C from 2MASS d Source for Vmag (0: PPMX (Tycho-2), 1: Transformed from USNO- 
B1.0, 2: Transformed from 2MASS) e Blending Flag (0: unblended, 1: unlikely blend, 2: possible blend, 3: probable 
blend, 4: problematic) 



TABLE 8 

Catalog of candidate periodic variable stars. II: AoV 





ID 


p a 

r ae 


SN ae a 


Amp aee b 


Amp aet b 


n a > c 


P * a 


SN at a 


Amp ate b 


Amp a ii b 


Q at a ' c 






d 




mag 


mag 




d 




mag 


mag 




HAT 


■085-0004278 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


HAT- 


■085-0053487 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


14.172864 


26.52 


0.412 


0.202 


1 


HAT- 


086-0001701 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


18.173471 


24.28 


0.010 


0.006 





HAT- 


086-0002357 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


HAT- 


086-0003825 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


HAT- 


■086-0004017 


54.656399 


79.24 


0.014 


0.004 


1 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


HAT- 


■086-0004248 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


HAT 


086-0005153 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


HAT- 


086-0005542 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


HAT- 


086-0005602 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


42.385695 


24.98 


0.015 


0.004 


1 



a The first 1 (or 2) letter(s) in the subscript denote the period-finding algorithm (a: AoV, ah: AoVHarm, d: DACF), the last letter 
in the subscript denotes the light-curve type (e: EPD, t: TFA). b The peak-to-peak amplitude. The first 2 (or 3) letter(s) in the 
subscript denote the period used, the last letter in the subscript denotes the light-curve type used to calculate the amplitude (e: 
EPD, t: TFA). C Quality flag (0: robust detection, 1: questionable detection, 2: not detected by the given method for the given 
light curve type). 



TABLE 9 

Catalog of candidate periodic variable stars. Ill: AoVHarm 





ID 


P u a 

r arte 


SN ahe a 


Amp ahee b 


Amp n i,et b 




Paht" 


SN aht a 


Amp ahte b 


Amp aM i b 


Qaht c 






d 




mag 


mag 




d 




mag 


mag 




HAT 


■085-0004278 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


HAT- 


■085-0053487 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


HAT- 


086-0001701 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


17.887814 


41.17 


0.010 


0.005 





HAT- 


086-0002357 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


HAT 


086-0003825 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


HAT- 


■086-0004017 


40.868328 


95.12 


0.014 


0.011 


1 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


HAT- 


■086-0004248 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


HAT 


086-0005153 


16.368924 


49.27 


0.010 


0.003 





-99.999999 


-999.99 


-9.999 


-9.999 


2 


HAT- 


086-0005542 


20.363047 


64.94 


0.010 


0.009 





20.411745 


29.90 


0.010 


0.009 





HAT- 


086-0005602 


-99.999999 


-999.99 


-9.999 


-9.999 


2 


21.071980 


47.14 


0.013 


0.005 


1 



a The first 1 (or 2) letter(s) in the subscript denote the period-finding algorithm (a: AoV, ah: AoVHarm, d: DACF), the last letter 
in the subscript denotes the light-curve type (e: EPD, t: TFA). b The peak-to-peak amplitude. The first 2 (or 3) letter(s) in the 
subscript denote the period used, the last letter in the subscript denotes the light-curve type used to calculate the amplitude (e: EPD, 
t: TFA). C Quality flag (0: robust detection, 1: questionable detection, 2: not detected by the given method for the given light curve 
type). 



K/M Dwarf Variability Survey 



TABLE 10 

Catalog of candidate periodic variable stars. IV: DACF EPD 





ID 


Pde a 


eP de a ' b 




Pr de a ' d 


Amp dee ° 


Amp det c 


Pfde a ' f 


Q de a ' s 








d 


d 






mag 


mag 






HAT- 


■085- 


■0004278 


28.878214 


0.744375 


0.2040 


-9.9122 


0.021 


0.006 





1 


HAT- 


■085- 


■0053487 


-99.999999 


-99.999999 


-9.9999 


999.9999 


-9.999 


-9.999 


2 


2 


HAT- 


-086- 


-0001701 


-99.999999 


-99.999999 


-9.9999 


999.9999 


-9.999 


-9.999 


2 


2 


HAT- 


■086- 


-0002357 


28.562554 


1.731345 


0.0389 


-21.1846 


0.009 


0.002 





1 


HAT- 


■086- 


■0003825 


28.930191 


0.703257 


0.1421 


-11.8251 


0.009 


0.002 


1 





HAT- 


-086- 


■0004017 


40.432676 


1.834407 


0.0633 


-12.3304 


0.013 


0.014 


1 


1 


HAT- 


■086- 


-0004248 


30.140252 


0.840410 


0.1182 


-12.9545 


0.011 


0.002 


1 


1 


HAT- 


■086- 


-0005153 


50.034209 


3.677269 


0.1127 


-13.5771 


0.007 


0.008 


1 





HAT- 


■086- 


-0005542 


-99.999999 


-99.999999 


-9.9999 


999.9999 


-9.999 


-9.999 


2 


2 


HAT- 


■086- 


-0005602 


64.435201 


11.503159 


0.0130 


-11.8146 


0.009 


0.005 


1 






a The first 1 (or 2) letter(s) in the subscript denote the period-finding algorithm (a: AoV, ah: AoVHarm, d: 
DACF), the last letter in the subscript denotes the light-curve type (e: EPD, t: TFA). b Period uncertainty 
from fitting a quadratic function to the cross-correlation peak. c Peak cross-correlation value d Natural loga- 
rithm of the false-alarm probability The peak-to-pcak amplitude. The first 2 (or 3) letter(s) in the subscript 
denote the period used, the last letter in the subscript denotes the light-curve type used to calculate the 
amplitude (e: EPD, t: TEA). Period flag (0: ok, 1: true period is likely a harmonic of measured period, 
2: no detection). 8 Quality flag (0: robust detection, 1: questionable detection, 2: not detected by the given 
method for the given light curve type). 



TABLE 11 

Catalog of candidate periodic variable stars. V: DACF TFA 



ID 


P di a 


eP dt a ' b 


y<« a ' c 


p rdt a,d 


Amp dte ° 


Amp dtt e 


Pfdt a ' f 


Q di a ' g 




d 


d 






mag 


mag 






HAT-085-0004278 


-99.999999 


-99.999999 


-9.9999 


999.9999 


-9.999 


-9.999 


2 


2 


HAT-085-0053487 


-99.999999 


-99.999999 


-9.9999 


999.9999 


-9.999 


-9.999 


2 


2 


HAT-086-0001701 


-99.999999 


-99.999999 


-9.9999 


999.9999 


-9.999 


-9.999 


2 


2 


HAT-086-0002357 


-99.999999 


-99.999999 


-9.9999 


999.9999 


-9.999 


-9.999 


2 


2 


HAT-086-0003825 


-99.999999 


-99.999999 


-9.9999 


999.9999 


-9.999 


-9.999 


2 


2 


HAT-086-0004017 


-99.999999 


-99.999999 


-9.9999 


999.9999 


-9.999 


-9.999 


2 


2 


HAT-086-0004248 


62.607406 


36.891590 


0.0150 


-10.7230 


0.014 


0.004 


1 


1 


HAT-086-0005153 


-99.999999 


-99.999999 


-9.9999 


999.9999 


-9.999 


-9.999 


2 


2 


HAT-086-0005542 


-99.999999 


-99.999999 


-9.9999 


999.9999 


-9.999 


-9.999 


2 


2 


HAT-086-0005602 


20.668377 


2.433455 


0.0079 


-16.9868 


0.013 


0.005 









a The first 1 (or 2) letter(s) in the subscript denote the period-finding algorithm (a: AoV, ah: AoVHarm, d: 
DACF), the last letter in the subscript denotes the light-curve type (e: EPD, t: TFA). b Period uncertainty 
from fitting a quadratic function to the cross-correlation peak. c Peak cross-correlation value d Natural loga- 
rithm of the false-alarm probability 6 The peak-to-peak amplitude. The first 2 (or 3) letter(s) in the subscript 
denote the period used, the last letter in the subscript denotes the light-curve type used to calculate the 
amplitude (e: EPD, t: TFA). f Period flag (0: ok, 1: true period is likely a harmonic of measured period, 
2: no detection). g Quality flag (0: robust detection, 1: questionable detection, 2: not detected by the given 
method for the given light curve type). 



Hartman et al. 



TABLE 12 

Catalog of candidate eclipsing binaries. I: Astrometric and Photometric Data 





ID 


RA a 


DE a 


pmRA b 


pmDE b 


Jmag c 


Hmag c 


Kmag c 


Vmag 


Vref 1 


Bl c 






h:m:s 


deg:m:s 


mas/yr 


mas/yr 


mag 


mag 


mag 


mag 






HAT- 


086-0022206 


23:50:17.08 


+51:11:28.9 


37.45 


-20.48 


11.655 


11.009 


10.858 


13.656 


1 


1 


HAT- 


087-0016301 


00:34:50.23 


+53:42:12.1 


-37.41 


12.76 


11.310 


10.613 


10.433 


13.904 


2 


2 


HAT- 


088-0017192 


01:38:56.50 


+48:52:40.2 


-11.17 


-29.81 


11.469 


10.807 


10.608 


13.920 


1 





HAT- 


113-0005699 


18:09:47.63 


+49:02:55.0 


-53.66 


-21.85 


11.362 


10.710 


10.560 


14.126 


1 





HAT- 


127-0008153 


03:04:05.20 


+42:03:10.5 


53.86 


-54.80 


10.650 


9.983 


9.730 


14.159 


1 


2 


HAT- 


127-0009895 


03:02:22.68 


+43:40:48.0 


-76.01 


-30.99 


11.018 


10.375 


10.198 


13.890 


1 


3 


HAT- 


129-0025342 


04:13:06.81 


+41:21:01.4 


42.67 


-26.33 


11.938 


11.372 


11.160 


14.311 


1 


1 


HAT- 


131-0026711 


05:16:36.90 


+48:35:44.3 


-14.58 


-32.73 


11.802 


11.138 


10.982 


14.275 


1 


2 


HAT- 


132-0012475 


05:43:20.18 


+47:17:00.1 


7.93 


-29.15 


11.345 


10.704 


10.569 


13.765 


1 


2 


HAT- 


133-0001901 


06:31:57.29 


+41:33:57.3 


-56.65 


-45.10 


9.581 


8.930 


8.787 


11.638 


1 


1 



a J2000, from 2MASS b from PPMX C from 2MASS d Source for Vmag (0: PPMX (Tycho-2), 1: Transformed from USNO- 
B1.0, 2: Transformed from 2MASS) e Blending Flag (0: unblended, 1: unlikely blend, 2: possible blend, 3: probable 
blend, 4: problematic) 



TABLE 13 

Catalog of candidate eclipsing binaries. II: Ephemerides 



ID 



Per 
d 



ePer 
d 



HJD0p a 
d 



eHJD0p b 
d 



DURp c 
h 



HJD0s a 
d 



eHJD0s b DURs c 
d h 



HAT 
HAT 
HAT 
HAT 
HAT 
HAT 
HAT 
HAT 
HAT 
HAT 



086- 0022206 

087- 0016301 

088- 0017192 
113-0005699 
127-0008153 
127-0009895 
129-0025342 

131- 0026711 

132- 0012475 

133- 0001901 



0.4290339 
3.1325961 
1.4030920 
0.2278737 
2.5681924 
0.5916726 
4.6109145 
0.6639531 
0.8663223 
0.2929328 



0.0000062 
0.0003763 
0.0000252 
0.0000013 
0.0000869 
0.0000239 
0.0002340 
0.0000048 
0.0000035 
0.0000001 



2454431. 
2454427. 
2453693. 
2453901. 
2452979. 
2452985. 
2453726. 
2454497. 
2454128. 
2453390. 



11547 
60691 
67973 
88979 
53695 
85734 
54374 
17095 
11503 
29599 



0.00042 
0.00504 
0.00066 
0.00028 
0.00140 
0.00187 
0.00129 
0.00026 
0.00145 
0.00016 



1.53 
5.24 
2.15 
1.24 
1.98 
4.16 
2.85 
1.43 
3.12 
0.91 



2454436. 
2454416. 
2453680. 
2453905. 
2452991. 
2452990. 
2453724. 
2454498. 
9999999. 
2453376. 



90565 
64515 
35434 
42128 
10638 
88691 
34418 
83198 
99999 
38208 



0.00117 
0.00442 
0.00071 
0.00039 
0.00215 
0.05334 
0.00235 
0.00170 
9.99999 
0.00022 



1.52 
4.63 
2.00 
1.14 
1.53 
4.55 
3.34 
1.42 
9.99 
0.88 



a Epoch of primary (p) or secondary (s) eclipse. The value is 9999999.99999 if a secondary eclipse is not detected. 15 Un- 
certainty on the epoch of eclipse. c Eclipse duration. 



TABLE 14 

Catalog of flare events. I: Astrometric and Photometric Data 3 - 





ID 


RA b 


DE b 


pmRA c 


pmDE c 


Jmag d 


Hmag d 


Kmag d 


Vmag 


Vrcf° 






h:m:s 


deg:m:s 


mas/yr 


mas/yr 


mag 


mag 


mag 


mag 




HAT- 


087-0002617 


00:21:57.81 


+49:12:38.0 


205.00 


-32.56 


9.139 


8.454 


8.205 


12.474 


1 


HAT 


■094-0007885 


05:35:22.84 


+53:16:02.7 


4.39 


-35.34 


10.797 


10.139 


9.958 


13.797 


1 


HAT 


101-0000851 


09:59:39.19 


+48:47:54.6 


-45.67 


8.38 


9.686 


9.002 


8.817 


12.137 


1 


HAT 


103-0001602 


11:03:17.04 


+49:48:25.5 


-63.32 


-7.74 


10.898 


10.286 


10.114 


13.544 


1 


HAT 


119-0031226 


22:02:08.40 


+49:07:17.3 


-68.41 


0.39 


11.102 


10.434 


10.228 


12.777 


1 


HAT 


123-0002064 


00:53:06.48 


+48:29:38.5 


220.82 


-162.32 


9.427 


8.773 


8.580 


12.899 





HAT 


123-0002064 


00:53:06.48 


+48:29:38.5 


220.82 


-162.32 


9.427 


8.773 


8.580 


12.899 





HAT 


123-0003289 


00:41:41.41 


+44:10:53.1 


-46.23 


-15.14 


10.012 


9.399 


9.188 


14.251 


1 


HAT 


126-0002065 


02:40:52.51 


+44:52:36.5 


215.74 


3.85 


9.277 


8.696 


8.461 


13.197 





HAT 


128-0007422 


03:42:45.17 


+43:27:40.1 


30.96 


-64.10 


10.435 


9.830 


9.568 


14.611 


1 



a Stars may be listed more than once if they have multiple flare events b J2000, from 2MASS C from PPMX d from 
2MASS 6 Source for Vmag (0: PPMX (Tycho-2), 1: Transformed from USNO-B1.0, 2: Transformed from 2MASS) 



K/M Dwarf Variability Survey 



TABLE 15 

Catalog of flare events. II: Flare data 8. 





ID 


HJD b 


Amp c 


r d 


Nobs 6 






d 


d 




HAT- 


087-0002617 


2454437.61429 


0.130 


0.0047 


26 


HAT- 


■094-0007885 


2454463.88452 


0.219 


0.0092 


7 


HAT- 


101-0000851 


2454162.63179 


0.200 


0.0176 


26 


HAT- 


103-0001602 


2454563.96193 


0.279 


0.0212 


29 


HAT- 


119-0031226 


2453360.74326 


0.552 


0.0346 


31 


HAT- 


123-0002064 


2454337.05672 


0.118 


0.0053 


11 


HAT- 


123-0002064 


2454391.64587 


0.482 


0.0065 


15 


HAT- 


123-0003289 


2454338.90758 


0.236 


0.0087 


23 


HAT- 


126-0002065 


2452999.82533 


0.277 


0.0056 


9 


HAT- 


128-0007422 


2453387.68872 


0.346 


0.0216 


13 



a Stars may be listed more than once if they have multiple 
flare events' 3 Time of flare peak c Peak intensity relative to the 
non-flaring intensity. d 1/e decay time of the flare. c Number of 
observations included in fitting the flare. 



TABLE 16 

Cross-matches to previously identified variable stars 



ID 


Alternate ID 


Mflag a 


BlFlag b 


HAT-086-0005602 


V0544 Cas 


1 


3 


HAT-086-0022206 


V1001 Cas 





1 


HAT-091-0010976 


V0524 Per 





2 


HAT-091-0016338 


V0485 Per 





1 


HAT-112-0001188 


1SWASP J173659.28+485946.1 








HAT-113-0005699 


ROTSE1 J180947.50+490254.5 








HAT-113-0005699 


VI 104 Her 








HAT-121-0003519 


NSV 14597 





3 


HAT-128-0006973 


NSV 01380 


1 


3 


HAT-128-0015056 


NSVS 15796 


1 


3 



a 0: the match is correct, 1: the two stars are not the same, but may be 
blended in the HAT images. b Automatically generated blending flag. 0: 
unblended, 1: unlikely blend, 2: potential blend, 3: probable blend, 4: 
problematic, 5: source is a flare star, no blending info is generated. 



TABLE 17 

Cross-matches to ROSAT X-ray sources 



ID 



RosatID a 



lgfxfjt 



HAT-087-0001873 1RXS J002854. 1+502226 -1.7955 

HAT-087-0002547 1RXS J003633. 7+553728 -2.2528 

HAT-087-0002617 1RXS J002158.0+491245 -2.1767 

HAT-093-0001111 1RXS J044853.2+552725 -2.3278 

HAT-094-0004894 1RXS J050718.8+530713 -2.3805 

HAT-094-0005389 1RXS J053001. 7+484932 -2.0147 

HAT-094-0006428 1RXS J050809.3+532437 -1.7462 

HAT-094-0009460 1RXS J050012.4+532033 -1.6652 

HAT-094-0384094 1RXS J053625. 1+533618 -2.2575 

HAT-103-0001602 1RXS J110318. 0+494822 -1.8031 



a ID of matching ROSAT source 13 Estimated log(/^//, 7 ) 
(ratio of X-ray flux to J-band flux) . 



